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EXECUTIVE  SUMMARY 


Program  SOM-LA  (Seat/Occupant  Model  -  Light  Aircraft)  has  been 
developed  for  use  in  evaluating  the  crashworthiness  of  aircraft 
seats  and  restraint  systems.  It  combines  a  three-dimensional  dy¬ 
namic  model  of  the  human  body  with  a  finite  element  model  of  the 
seat  structure.  It  is  intended  to  provide  the  design  engineer  a 
tool  with  which  he  can  analyze  the  structural  elements  of  the  seat 
as  well  as  evaluate  the  dynamic  response  of  the  occupant  durinq  a 
crash. 

The  occupant  model  consists  of  twelve  masses  that  represent 
the  upper  and  lower  torso,  neck/  head/  and  two  segments  for  each 
of  the  arms  and  legs.  An  optional  model  of  the  human  body  includes 
beam  elements  in  the  spine  and  neck,  but  is  restricted  to  two- 
dimensional  motion. 

External  forces  are  applied  to  the  occupant  by  the  cushions, 
floor,  and  restraint  system.  Interface  between  the  occupant  and 
seat  is  provided  by  the  seat  bottom  cushion,  back  cushion,  and  an 
optional  headrest.  The  restraint  system  can  consist  of  a  lap  belt 
alone  or  combined  with  a  single  shoulder  belt,  over  either  shoul¬ 
der,  or  a  double-strap  shoulder  harness.  A  lap  belt  tiedown 
strap,  or  negative-G  strap,  can  also  be  included.  Each  component 
of  the  restraint  system  can  be  attached  to  either  the  seat  or  the 
aircraft  structure. 

For  two  standard  occupants,  a  5 Oth-percentile  human  male  and 
a  50th-percentile  (Part  572)  anthropomorphic  dummy,  all  dimensions 
and  inertial  properties  required  for  simulation  are  stored  within 
the  program.  For  other  occupants,  these  characteristics  must  be 
provided  as  input. 

The  seat  structure  is  modeled  using  the  finite  element  method 
of  analysis,  selected  because  it  is  not  dependent  on  previous  test¬ 
ing,  and  it  has  the  flexibility  to  deal  with  a  wide  range  of  de¬ 
sign  concepts.  The  SOM-LA  seat  analysis  includes  triangular  plate 
elements,  three-dimensional  beam  elements,  and  spring  elements. 

It  has  the  capability  to  model  large  displacements,  nonlinear  ma¬ 
terial  behavior,  local  buckling,  and  various  internal  releases  for 
beam  elements. 

The  digital  computer  program  based  on  the  occupant  and  seat 
models  described  above  has  been  written  entirely  in  FORTRAN  to  en¬ 
sure  a  high  degree  of  compatibility  with  various  digital  computer 
systems.  During  development,  the  program  has  been  run  on  IBM, 
Univac,  and  CDC  computer  systems.  Output  data  include  occupant 
segment  positions,  velocities,  and  accelerations;  restraint  system 
and  cushion  forces;  injury  criteria;  and  details  of  contact  be¬ 
tween  the  occupant  and  the  aircraft  interior.  Seat  output  includes 
nodal  displacements,  element  stresses,  and  forces  at  the  points 
of  attachment  to  the  aircraft  structure. 
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Validation  has  been  based  on  data  from  several  series  of  de¬ 
celeration  sled  tests  conducted  at  the  FAA  Civil  Aerom-'dical  In¬ 
stitute  (CAMI).  The  response  of  the  combined  occupant  and  seat 
models  was  verified  by  comparison  with  data  from  tests  that  uti¬ 
lized  specially  designed  and  fabricated  seats  with  replaceable 
legs.  The  test  conditions  in  that  series  were  specifically  se¬ 
lected  to  cause  significant  plastic  deformation  of  the  legs.  Re¬ 
sponse  of  the  occupant  model,  particularly  to  a  vertical  input 
acceleration,  was  validated  using  data  from  other  test  series 
conducted  with  a  rigid  seat  and  with  a  production  energy-absorbing 
helicopter  seat.  The  final  phase  of  validation  included  simula¬ 
tion  of  dynamic  tests  of  production  general  aviation  seats.  The 
tests  are  described  in  the  report,  and  model  predictions  are  com¬ 
pared  with  test  data.  The  SOM-LA  program  is  shown  to  be  a  poten¬ 
tially  useful  tool  for  analyzing  the  performance  of  a  seat  and  re¬ 
straint  system  in  a  crash  environment. 


1.0  INTRODUCTION 


The  design  of  crashworthy  seats  and  restraint  systems  for  aircraft 
presents  a  complex  engineering  problem,  the  solution  of  which  can 
be  greatly  aided  by  sufficiently  rigorous  analytical  techniques. 

The  crash  environment  can  vary  widely  from  one  accident  to  another, 
thus  a  great  number  of  conditions  must  be  evaluated  to  establish 
those  critical  to  occupant  survival.  For  example,  the  restraint 
system  must  limit  the  movement  of  the  occupant  sufficiently  to 
eliminate  the  possibility  of  head  strike  on  rigid  cockpit  struc¬ 
ture.  Also,  the  relatively  low  tolerance  of  the  human  body  to 
accelerations  in  a  direction  parallel  to  the  spine  requires  the 
consideration  of  vertical  impact  forces  which  are  usually  present 
and  often  significant  in  crashes  of  light  fixed-wing  aircraft  and 
helicopters  (reference  1).  A  very  strong,  rigid  seat  is  not  a 
truly  valid  solution,  since  it  would  not  only  incur  serious  weight 
penalties,  but  would  transmit  high  vertical  impact  forces  directly 
to  the  occupant.  In  light  aircraft  and  helicopters,  it  is  seldom 
practical  to  consider  designing  sufficient  energy-absorbing  capabil¬ 
ity  into  the  lower  airframe  structure  to  protect  against  these 
vertical  forces,  as  the  crush  space  is  generally  not  available. 
Rather,  a  crashworthy  seat  for  these  aircraft  should  include  the 
capacity  to  absorb  energy  through  controlled  deformation  in  the 
vertical  direction,  thus  reducing  the  accompanying  loads. 

In  the  initial  design  phases,  it  is  desirable  to  evaluate,  in  some 
detail,  existing  seats  and  restraint  systems  in  their  surrounding 
cockpits,  thus  establishing  existing  weaknesses.  It  is  then  de¬ 
sirable  to  make  modifications  and  to  evaluate  the  effect  of  these 
modifications  on  improving  the  survivability  of  the  system.  These 
evaluations  must  be  conducted  for  a  great  many  of  the  possible 
crash  environments,  thus  constituting  a  relatively  large  matrix. 
Testing  is  extremely  expensive  and  requires  a  great  deal  of  time, 
since  design  modifications  must  be  developed  and  fabricated  prior 
to  testing.  Therefore,  an  analytical  technique,  such  as  was  devel¬ 
oped  in  this  program,  is  required. 
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A  number  of  one-,  two-,  and  three-dimensional  mathematical  models 
of  the  human  body  have  been  developed  for  crash  survivability 
analysis.  These  models  vary  in  complexity  and  possess  from  one  to 
forty  degrees  oi  freedom.  The  simplest  models  have  been  devel¬ 
oped  primarily  for  prediction  of  injury  to  a  single  component  or 
subsystem  of  the  body,  such  as  vertebral  injury  due  to  the  verti¬ 
cal  force  that  might  be  experienced  in  a  crash  or  the  firing  of 
an  ejection  seat.  Simulation  of  the  three-dimensional  response 
of  the  entire  body  requires  many  more  degrees  of  freedom,  but  per¬ 
mits  more  general  use.  Most  of  the  three-dimensional  models  have 
been  developed  for  use  in  evaluation  of  automobile  interior  de¬ 
sign  with  respect  to  injuries  caused  by  secondary  impacts,  such 
as  the  three-dimensional  models  described  in  references  2  through 
5.  Seats  have  been  represented  in  a  very  simple  manner  because 
i.i  automobiles  the  role  of  the  seat  design  in  determining  occu¬ 
pant  survival  is  minimal.  Therefore,  a  simulation  model  intended 
specifically  for  the  aircraft  application  is  required. 

The  development  oc  a  chree-dimensional  mathematical  model  of  a 
light  aircraft  seat,  occupant,  and  restraint  system  is  described 
in  this  report.  This  model  forms  the  basis  for  a  simulation  com¬ 
puter  program  that  has  been  written  specifically  for  use  in  crash- 
worthy  design  and  analysis  of  light  aircraft  seats  and  restraint 
systems,  Program  SOM-LA  (Sont/Cccupant  Model-Light  Aircraft).  The 
program,  which  combines  a  finite  olement  model  of  the  soat  struc¬ 
ture  with  a  lumped-parameter  model  of  the  aircraft  occupant,  has 
boon  organized  so  as  to  minimise  the  volume  and  complexity  of  in¬ 
put  data  and  to  focus  on  seat  and  rertraint  system  design  param¬ 
eters. 


This  original  mcdel  was  described  in  a  comprehensive  technical 
report  that  was  published  by  the  FAA  m  1975  {reference  6).  A 
number  of  modifications  have  been  made  to  the  medei  since  then  to 
improvo  simulation  quality  and  to  prov '.ae  increased  capability 
and  additional  desirable  output.  In  1977,  a  testing  program  was 
initiated  by  the  Civil  Aeroraedical  Institute  (CAMI)  to  provide 
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data  for  validation  of  the  model.  The  validation  effort  pointed 
out  areas,  particularly  in  the  seat  structure  model,  where  modi¬ 
fications  to  the  program  were  required.  Also,  use  of  the  program 
both  by  its  developers  and  other  organizations  demonstrated  the 
need  for  additional  capabilities.  Probably  the  most  significant 
change  to  the  program  since  publication  of  reference  6  has  been 
the  incorporation  of  a  more  general  finite  element  seat  model  with 
capabilities  of  simulating  large,  inelastic  deformations  and  crip¬ 
pling  of  hollow  members.  With  respect  to  the  occupant  model,  the 
program  user  now  has  the  option  of  either  a  three-dimensional  oc¬ 
cupant  made  up  of  rigid  segments  or  a  two-dimensional  model  which 
has  beam  elements  representing  the  spine  and  neck,  thus  providing 
greater  capability  for  injury  prediction. 

Presented  in  this  report  are  the  details  of  the  occupant  and  seat 
models,  as  well  as  the  results  of  model  validation.  Operation  of 
the  computer  program  is  described  in  Volume  II  of  this  report. 
Program  SQM-LA  User  Manual. 


2.0  OCCUPANT  MODEL 


The  mathematical  model  includes  a  lumped-parameter  representation 
of  the  vehicle  occupant  and  a  finite  element  seat.  Interface  be¬ 
tween  the  seat  and  occupant  is  provided  by  seat  cushions  and  a  re¬ 
straint  system,  which  consists  of  a  lap  belt  and,  if  desired,  a 
single-strap  or  double-strap  shoulder  harness.  A  lap  belt  tiedown 
strap,  or  negative-G  strap,  can  be  added  for  simulation  of  a  full, 
five-point  restraint  system.  The  response  of  the  occupant  and 
seat  can  be  predicted  for  any  given  set  of  aircraft  impact  con¬ 
ditions,  including  the  initial  velocity  and  attitude  and  the  in¬ 
put  acceleration. 

As  mentioned  in  chapter  1.0,  both  a  three-dimensional  occupant 
model  consisting  of  rigid  links  and  a  two-dimensional  model  with 
deformable  spinal  elements  are  included.  This  chapter  provides  a 
discussion  of  the  development  of  both  models,  including  details 
of  the  approach  to  formulating  the  equations  of  motion  and  of 
the  technique  used  for  their  solution.  Section  2.1  presents  the 
three-dimensional  model  and  section  2.2,  the  two-dimensional  model. 
Subsequent  sections  cover  aspects  of  the  equations  which  apply  to 
both  models,  specifically  the  body  joint  model,  the  treatment  of 
external  forces,  and  body  dimensions  and  properties. 

2 . 1  THREE-DIMENSIONAL  OCCUPANT  MODEL 


The  three-dimensional  mathematical  model  of  the  aircraft  occupant 
is  made  up  of  twelve  rigid  segments,  as  shown  in  figure  1.  This 
number  is  thought  to  represent  the  minimum  that  will  permit  ac¬ 
curate,  meaningful  simulation  of  three-dimensional  response.  A 
greater  number  might  possibly  improve  the  accuracy  of  simulation 
but  would,  in  turn,  increase  program  execution  cost.  Arm  and  leg 
segments  are  included  to  enable  prediction  of  injuries  to  these 
extremities.  Although  leg  and  arm  injuries,  in  themselves,  may 
not  be  as  serious  as  head  or  chest  injuries,  they  may  prevent  ex¬ 
cape  from  a  stricken  aircraft  and  the  potential  hazard  of  post¬ 
crash  fire. 
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$  Segment  mass  center 


Figure  1.  Twelve -segment  (three-dimensional) 
occupant  model. 

Each  of  the  body  joints,  with  the  exception  of  the  elbow,  knee, 
and  head-neck  joints,  possesses  three  rotational  degrees  of  free¬ 
dom.  Because  of  the  hinge-type  motion  of  a  forearm  or  lower  leg 
relative  to  an  upper  arm  or  thigh,  respectively,  the  position  of 
each  of  these  segments  is  described  by  one  additional  angular  co¬ 
ordinate,  as  is  the  position  of  the  head  relative  to  the  neck. 
Therefore,  the  occupant  system  possesses  a  total  of  29  degrees  of 
freedom. 


2.1.1  Coordinate  Systems 


Fixed  at  the  center  of  mass  of  each  of  the  12  segments  is  a 
right-handed  Cartesian  coordinate  system.  For  segment  n  (n  =  1 
2,...,  12)  the  local  coordinate  system  is  denoted  by  axes  (xn, 
y  ,  2  ) .  Positive  directions  are  defined  such  that  when  the  body 
is  seated  as  shown  in  figure  2,  with  the  torso  and  head  upright, 
the  upper  arms  parallel  to  the  torso,  and  the  elbows  and  knees 
bent  at  right  angles,  positive  x^  is  directed  forward,  yR  to  the 
left,  and  zn  upward. 

In  order  to  describe  a  general  position  of  the  body,  it  is  nec¬ 
essary  to  relate  the  orientation  of  each  segment  (xn,  yn,  zn)  to 
the  inertial  system  (X,  V,  Z).  The  angular  relationship  between 
the  local,  segment-fixed  coordinates  and  the  inertial  system  can 
be  expressed  by  the  transformation 


Because  three  angular  coordinates  can  be  used  to  define  the  rota¬ 
tion  of  a  given  segment,  it  is  convenient  to  utilize  a  set  of  co¬ 
ordinates  that  will  suffice  as  generalized  coordinates  in  the  for¬ 
mulation  of  tho  equations  of  motion.  A  system  of  Eulerian  angles 
provides  a  convenient  set  of  three  independent  angular  coordinates. 
Assuming  that  the  local  (xn,  yn,  zn)  system  is  initially  coinci¬ 
dent  with  the  inertial  (X,  Y,  Z)  system,  the  Euler  angles  are  a 
scries  of  three  rotations,  which,  when  performed  in  the  proper  se¬ 
quence,  permit  the  system  to  attain  any  orientation  and  uniquely 
define  that  position.  The  particular  set  of  Euler  angles  selected 
for  use  here  is  illustrated  in  figure  3  and  defined  as  follows: 


Figure  3.  The  Euler  angles. 


1.  A  positive  rotation  i|»  about  the  Z-axis,  resulting  in 
the  primed  (x',  y' ,  z1)  system. 

2.  A  positive  rotation  0  about  the  y'-axis  resulting  in 
the  double-primed  (x*,  y",  z")  system. 

3.  A  positive  rotation  ^  about  the  x*-axis  resulting  in 
the  final  (x#  y,  z)  system. 

In  order  to  determine  the  elements  of  the  transformation  matrix 
(Tn],  it  is  necessary  to  consider  the  matrix  equations  that  in¬ 
dicate  the  three  individual  rotations  previously  described  above. 
Referring  again  to  these  definitions  of  ^ ,  9,  and  $ ,  the  follow¬ 
ing  equations  are  obtained: 
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{r*  }  «  td>)  {r ) 
or 

(R}  *  U)  (B)  U)  {r)  (5) 

whore  (R!  represents  the  components  of  a  vector  in  the  inertial 
system  and  {rl  represents  the  same  vector  in  the  final  (x,  y#  z) 
system.  Performing  the  matrix  multiplications  indicated  in  equa¬ 
tion  ( 5 ) ,  the  elements  of  the  transformation  matrix  in  equation  (1) 
are  obtained: 
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The  additional  constraint  of  hinge-type  rotation,  at  the  elbows 
and  knees,  requires  the  use  of  one  additional  angular  coordinate 
to  define  the  position  of  each  of  the  forearm  and  lower  leg  seg¬ 
ments.  Referring  to  figure  4,  the  angular  position  of  the  fore¬ 
arm  segments  (  t  a  5,  7)  is  given  by 
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and  the  lower  leg  segments  (m  *  9,  11)  by 
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From  equations  (6)  and  (7),  the  elements  of  the  transformation 
matrix  for  the  forearms  (segments  S  and  7)  arc  written  as: 
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Figure  4. 


Definition  of  angular  coordinates 
a  for  elbows,  knees,  and  neck. 
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prow  equations  (6)  and  ( 8 ) ,  the  elements  of  the  transformation 
matrix  elements  for  the  logs  (segments  9  and  11)  are  obtained: 
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Segment  12,  the  neck,  was  added  later.  The  hinge-type  joint  be¬ 
tween  the  neck  and  the  head  is  treated  similarly  to  the  knee  joints, 
adding  the  generalized  coordinate  9^/  as  illustrated  in  figure  4. 

Having  developed  the  relationships  expressed  in  the  equations  (1) 
through  (10),  the  position  of  the  occupant  can  be  described  by  the 
following  set  of  generalized  coordinates: 
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The  above  coordinates  include  the  Cartesian  coordinates  of  the 
mass  center  of  segment  1  (X^,  Y^,  Z^),  selected  as  a  reference 
point  on  the  body,  seven  sets  of  Eulerian  angles,  and  the  five 
additional  angular  coordinates  for  the  elbows,  knees,  and  neck. 
Positions  of  the  segment  mass  centers  are  presented  in  appendix  A. 

2.1.2  Lagrange *s  Equations 

The  response  of  the  occupant  system  is  described  by  Lagrange's 
equations  of  motion,  which  are  written  for  the  29  generalized  co¬ 
ordinates.  The  equations  are  developed  according  to 

fr  <— )  -  “  “  <3  “  2,...,  29)  (12) 

aq  j  5qj  3 

whore  L  is  tho  Lagrangian  function 

L  «  T  -  V  (13) 

t  represent r  time,  are  the  generalized  forces  not  derivable 
from  a  potential  function.  (Forces  that  are  derivable  from  a  po¬ 
tential  function  are  obtained  from  L,  and  T  and  V  are  tho  system 
kinetic  and  potential  energies,  respectively.) 


14 


Because  the  system  being  treated  does  not  involve  any  velocity- 
dependent  potentials,  equation  (12)  can  be  written  as 


d3  , 3T  .  3T 
dt  a  *  3q  . 

3qj 


( j  B  1»  2, . . . ,  29) 


(14) 


The  system  kinetic  energy  contains  both  translational  and  rotational 
parts: 
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where  M  is  the  mass  of  segment  n  and  1,1,  and  I  are  mass 
n  x  y  2 

moments  of  inertia  of  segment  n  with  respect  to  the  local  coordi¬ 
nate  axes  (xn,  yn,  zn),  assumed  to  be  principal  moments  of  iner¬ 
tia. 


The  absolute  velocities  of  the  12  mass  segments  required  for  the 
translational  kinetic  energy  must,  of  course,  be  written  as  func¬ 
tions  of  the  generalized  coordinates  and  generalized  velocities 
in  order  to  use  equation  (14).  The  angular  velocity  components 

(«.>  *  ui  #  w  )  seen  in  equation  (15)  are  parallel  to  the  local 

n  yn  zn 

(x  ,  yn#  z  )  coordinate  systems.  These  angular  velocity  compo¬ 
nents  cannot  be  used  directly  in  Lagrange's  equations  because  they 
do  not  correspond  to  the  time  derivatives  of  any  set  of  coordinates 
that  specify  the  position  of  the  segment.  They  must  be  written  as 
functions  of  the  generalized  coordinates,  U9ing  the  generalized 
angular  velocities  (|n,  0  ,  !n),  which  are  parallel  to  the  axes  2, 
y^,  and  xf",  respectively. 
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An  arbitrary  angular  velocity  of  segment  n,  w  ,  can  be  expressed 
as  a  function  of  the  generalized  angular  velocities  according  to 


«n  =  +  2n  +  *n 


(16) 


Referring  to  figure  3,  $n,  e  ,  and  do  not,  in  general,  form  a 
mutually  perpendicular  vector  triad.  and  £  are  both  perpendi- 
cular  to  £.  but  are  not  necessarily  perpendicular  to  each  other.) 
However,  they  can  be  considered  as  a  nonorthogonal  set  of  compo¬ 
nents  of  w  since  their  vector  sum  is  equal  to  w.  Summing  the  or¬ 
thogonal  projections  of  ,  §  ,  and  |n  on  the  (xR,  yn,  zR)  axes 
yields  the  angular  velocity  exponents  required  for  the  kinetic 
energy  expression; 


“x  =  -  K  sln  8n 
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W  *  cos  en  sin  <Jn  ♦  0n  cos  <j>n 


(17) 


u>  =  i  cos  6  cos  A  -  O  sin  <j> 
z  n  n  n  n  n 

n 

The  system  potential  energy  is  simply  gravitational  potential, 
which  is  written  as 


12 

V  -  r  M  g  (Z  -  ) 

nil  n  n  no 


(18) 


whore  g  is  the  acceleration  due  to  gravity  and  Z  is  an  arbitrary 

no 

datura, 

2.1*3  Matrix  Equations 

for  the  purposes  of  computation,  the  equations  of  motion  are  re¬ 
written  in  the  following  form; 


lA(q)]  iq)  *  (D(q,q)l  ♦  (P(q)l  ♦  (R(q,q))  ♦  (Q(q,q)l  (19) 
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where  the  elements  of  the  inertia  matrix  [A]  and  the  vector  { B} 
are  derived  from  the  kinetic  energy  derivatives  of  Lagrange's 
equations.  In  other  words. 


*  *  • 

“  Bj  (q^,  q2 » •  •  •  /  <?29/  '  ^2 '  *  *  * '  ^29  ^ 

(j  ~  1,  2,...,  29)  (20) 

The  force  vector  (p)  is  derived  from  the  system  potential  energy 
according  to 

3  V 

Pj  (q  1 ,  q2>**,»  ^29*  s  "*  gq  (3  “  2,...,  29)  (21) 

J 

Both  { R )  and  (Q)  are  vectors  of  generalized  forces  derived  from  the 
right-hand  side  of  Lagrange's  equations.  The  vector  (r)  describes 
the  resistance  of  the  body  joints  to  rotation,  discussed  in  detail 
in  section  2.3;  (Q)  is  the  vector  of  generalized  external  forces, 
discussed  in  detail  in  section  2*4. 

2.2  TWO-DIMENSIONAL  OCCUPANT  MODEL 


In  order  to  achieve  economical  program  solutions  for  cases  where 
occupant  response  is  expected  to  be  symmetrical  with  respect  to 
the  X-z  plane,  a  two-dimensional  occupant  model  option  was  in¬ 
cluded  in  Program  SOM-LA.  This  two-dimensional  model  vas  con¬ 
figured  like  the  three-dimensional  model  shown  in  figure  1,  with 
the  exception  of  all  joints  being  hinge-type  joints.  Because  of 
the  potential  for  vertebral  injury  in  aircraft  accidents  that  in¬ 
volve  a  significant  vertical  component  of  impact  velocity,  some 
measure  of  vertebral  loading  was  considered  desirable  in  the  oc¬ 
cupant  model.  The  two-dimensional  occupant  model  was  configured 
to  include  beam  elements  in  both  the  torso  and  neck,  as  shown  in 
figure  5,  replacing  joints  that  exist  in  the  three-dimensional 
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Figure  5.  Eleven-segment  (two-dimensional) 
occupant  model . 


model.  The  two-dimensional  model  has  11  degrees  of  freedom,  as 
illustrated  in  figure  6.  Simulation  of  purely  vertical  (♦G_)  im- 
pact  with  the  three-dimensional  occupant  model  produced  less  bend¬ 
ing  of  the  torso  than  observed  in  tests  with  either  dummies  or 
human  cadavers.  Therefore,  in  order  to  produce  the  moments  on 
the  vertebral  column  that  are  induced  by  +G2  acceleration  in  both 
human  and  dummy  occupants,  the  mass  centers  of  the  torso  segments. 
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segments  1  and  2,  are  offset  from  the  vertebral  column  by  the  di¬ 
mensions  e^  and  respectively.  The  dimension  e^,  by  which  the 
center  of  mass  of  the  head  segment  is  placed  forward  of  the  neck, 
is  used  in  both  three-dimensional  and  two-dimensional  models. 

Development  of  the  equations  of  motion  for  the  two-dimensional 
occupant  follows  a  procedure  similar  to  that  described  in  sec¬ 
tion  2.1  for  the  three-dimensional  model.  However,  the  procedure 
is  simplified  because  the  transformation  between  a  local,  segment- 
fixed  coordinate  system  and  the  inertial  system  is  dependent  only 
on  the  angle  8.  Transformations,  therefore,  take  the  form  of 
equation  (3).  The  position  of  the  center  of  mass  of  segment  2, 
with  respect  to  that  for  segment  1  and  the  position  of  the  center 
of  mass  of  segment  3,  relative  to  segment  2,  depend  on  the  length 
and  curvature  of  the  beam  elements  in  the  spine  and  neck,  respec¬ 
tively. 

A  detailed  representation  of  the  kinematics  and  mechanics  of  the 
spinal  column  would  be  complicated.  The  presence  of  the  interver¬ 
tebral  disks  leads  to  high  flexibility  in  bending  and  high  stiff¬ 
ness  in  compression.  In  addition,  the  column  is  tapered  and  pos¬ 
sesses  considerable  initial  curvature.  A  simplified  model  is  in¬ 
corporated  into  SOM-LA  by  using  continuous  beam  elements  for  the 
neck  and  spine.  These  beam  elements  are  intended  to  model  the 
flexural,  as  well  as  axial,  motion  of  the  spine  and  neck  and  are 
subject  to  the  following  assumptions.  The  deformed  beam  elements 
take  the  shape  of  circular  arc  segments,  therefore  assuming  flex¬ 
ure  to  be  primarily  due  to  the  applied  bending  moment.  This  is 
equivalent  to  the  assumption  in  beam  theory  that  the  span/depth 
ratio  is  large  and  that  the  modulus  of  elasticity  and  shear  mod¬ 
ulus  are  of  the  same  order  of  magnitude  so  that  shear  deformation 
is  negligible. 

As  illustrated  in  figure  7,  The  position  of  one  end  relative  to 
the  other  in  terms  of  arc  length  and  angular  coordinates  at  each 
end  of  the  arc,  is  then  given  by 
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Figure  7.  Position  of  segment  2  relative 
to  segment  1. 


Positions  of  the  mass  centers  for  the  eleven  body  segments  of  the 
two-dimensional  model  are  presented  in  appendix  B. 


2.3  JOINT  RESISTANCE 


The  form  of  the  joint  resistance  vector  {R>  in  equation  (19)  de¬ 
pends  on  the  user’s  selection  of  occupant  type  -  either  dummy  or 
human.  Although  both  joint  models  contain  the  same  types  of  el¬ 
ements,  a  nonlinear  torsional  spring  and  a  viscous  torsional  damp¬ 
er,  the  relative  contributions  of  each  of  these  elements  determine 
the  type  of  occupant. 

The  11  body  joints  for  the  three-dimensional  model,  illustrated 
in  figure  1,  are  defined  as  follows: 

joint  1  -  Back,  between  12th  thoracic  and  1st  lumbar 
vertebrae 

Joint  2  -  Torso-neck,  between  7th  cervical  and  1st  thoracic 
vertebrae 

Joint  3  -  Right  shoulder 
Joint  4  -  Right  elbow 
Joint  5  -  Left  shoulder 
Joint  6  -  Left  elbow 
Joint  7  -  Right  hip 
Joint  8  -  Right  knee 
Joint  9  -  Left  hip 
Joint  10  -  Left  knee 

Joint  11  -  Head-neck,  at  occipital  condyles 

The  angular  displacement  of  joint  i  from  its  reference  position 
(figure  2)  is  given  by  3..  If  (  i  ,  J  ,  km)  and  (i„  j',  k  )  are 

triads  of  unit  vectors  in  the  local  coordinate  systems  of  two  ad¬ 
jacent  segments  connected  at  joint  i,  as  shown  in  figure  8,  the 
joint  angle  is  given  by 

3i  »  cos"1  (km  •  kn)  (23) 
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Figure  8.  Joint  angle  between  segments  m  and  n. 


where  {  ^  .  kn)  is  the  scalar  product.  Considering  the  geometry 
of  the  occupant  model  in  the  reference  position,  the  8^  for  the 
11  joints  are  given  by 


=  cos  1  (k  ^  •  k 

ft2  »  cos"1  (i'2  *  ^12^ 
fl3  *  cos"1  (k2  *  kj) 


h 


«  cos 


-1 


k6> 
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I 


6  =  a 

P6  7 


=  cos"*1  (•t1  •  ig) 


B8  =  a9 


cos”1  U1  *  i1Q) 


B10  =  all 


=  cos”1  Ug  •  i:12) 


(24) 


If  at  each  joint  i,  a  moment  and  a  torsional  damper  with  co¬ 
efficient  act  to  resist  motion  of  the  joint,  then  the  virtual 
work  done  on  the  system  as  each  joint  i  undergoes  a  virtual  dis¬ 
placement  'S  is 


11 


aw  -  -  C  <Mi  £IV 

i»l 


(25) 


Since  tho  ft,  are  functions  of  the  generalized  coordinates  q^,  the 
virtual  displacements  Si^  can  be  expressed  in  terms  of  correspond* 
ing  virtual  displacements  of  the  q^.  In  general,  such  an  expres¬ 
sion  would  take  the  form 


AR  .  a 

1 


29  3ft. 

E  IT 

j.l  *lj  3 


(i-1,  2,...,  11) 


(26) 


where  tho  partial  derivatives  3 /  ^<3 ^  are  functions  of  the  gen¬ 
eralized  coordinates.  Substituting  into  equation  (2$)  gives 


aw 


i\W  a  - 


11  29  .  9ft, 

E  E  ("i  ♦  ->i  V  nr  S 

i*l  j-1  1  1  1  **3  1 


(27) 
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Changing  the  order  of  summation#  equation  (27)  can  be  written  in 
the  general  form 


29 

6W  =  £  R,  Sq.  (28) 

j=l  3  3 

where  R ^  are  the  generalized  forces  acting  on  the  system.  As 
seen  in  equation  (19),  the  generalized  forces  are  treated  as  two 
distinct  types:  joint  resistance  forces  and  external  forces. 
Since  the  joint  resistance  terms  are  being  treated  here,  the  gen¬ 
eralized  joint  forces  referred  to  as  R ^  will  be  considered  alone. 
Equation  (28)  becomes,  more  specifically 
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<*w  =  E  R,  6q  (29) 

j=l  3  3 


and,  from  equation  (27),  R^  can  be  written 


As  mentioned  earlier  in  this  section,  the  type  of  occupant  is  de¬ 
termined  by  the  relative  contributions  of  H4  and  J,  to  the  R, 
terms,  for  the  dummy  joint,  the  resisting  torque  is  constant 
throughout  the  normal  range  of  joint  motion  and  increases  rapidly 
along  a  third-order  curve  to  a  higher  value  at  the  limiting  dis¬ 
placement  s*  ,  as  shown  in  figure  9.  The  normal  valuos  are 

set  equal  to  those  resulting  from  the  joint-tightening  procedure 
of  SAK  Recommended  Practice,  Anthropomorphic  Test  Device  for 
Dynamic  Testing  -  SAB  *1963.  That  is,  the  body  joints  will  just 
support  a  1-G  load  in  the  reference  (seated)  position,  with  the 
exception  of  the  torso  joints,  which  will  support  a  2-G  load. 

In  addition  to  a  small  viscous  damping  term  with  constant 
Ji  is  included  for  energy  dissipation. 


Angular  displacement  -  0^ 


Figure  9.  Dummy  joint  resisting  torque. 


The  resistance  of  each  human  joint  consists  of  up  to  three  terms. 
The  primary  resisting  force  during  normal  joint  rotation  is  a  vis¬ 
cous  damping  term  with  constant  coefficient  In  a  manner  simi¬ 

lar  to  the  case  of  the  dummy ,  a  resisting  torque  is  applied  at  the 
limit  of  the  joint  range  of  motion,  as  shown  in  figure  10(a).  An 
additional  term  used  to  simulate  muscle  tone  is  the  moment  H*, 
which  drops  to  rero  after  a  small  angular  displacement  from  the 
initial  position,  provided  that  the  crash  deceleration  is  suffi¬ 
cient  to  overcome  it  (figure  10(b)). 

For  the  two-dimensional  model,  the  equations  presented  in  this 
section  are  correct  except  that  summations  are  performed  for  8 
joints  and  11  degraee  of  freedom,  rather  than  the  11  and  29  pre¬ 
sented  here,  respectively. 


(a) 


Displacement  from  initial 
position  -  (0,  -  ) 

1  *0 


(b) 


Figure  10.  Human  joint  resisting  torques:  { 
limiting  moment  and  (b)  muscular 
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CO 


a)  displacement- 
resistance. 
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2.4  EXTERNAL  FORCES 


The  vector  of  generalized  external  forces  {Q}  is  developed  in  a 
manner  similar  to  that  discussed  in  the  previous  section  for  the 
joint  resistance  vector.  Equations  presented  in  this  section  in¬ 
clude  summations  whose  upper  limits  are  correct  for  the  three- 
dimensional  model,  with  12  segments  and  29  degrees  of  freedom. 

The  approach  for  the  two-dimensional  model  is  the  same,  except 
for  the  use  of  11  segments  and  11  degrees  of  freedom. 

The  resultant  external  force  E\  acting  on  segment  i  is  given  by 


F,  =  Fx  i  +  F  j  +  F  k  (31) 

1  xi  yi  zi 


where  Pv  ,  Fv  ,  and  F_  are  components  in  the  inertial  (X,  Y,  Z) 
Xi  Yi  Zi 

system.  The  absolute  position  of  the  point  Pi  on  segment  i, 
where  the  resultant  force  acts,  can  be  represented  by 


rp  -  Xp  i  +  Yp  j  *  Zp  k  (32) 

pi  Pi  pi  Pi 


As  the  resultant  force  applied  to  each  segment  i  undergoes  a  vir¬ 
tual  displacement  Arp  ,  having  components  (§Xp  ,  SYp  ,  6Zp  ),  the 

i  i  i  i 

virtual  work  on  the  system  done  by  the  F^  is 
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AW  *  £  (Fy  £Xp  ♦  Fy  6 Yp  ♦  F»  6Zp  )  (33) 

i  =  l  xi  *i  i  i  Pi 


Writing  the  virtual  displacement  components  in  terms  of  the  ganer 
alized  coordinates  q^: 
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(34) 


The  contact  forces  applied  to  the  occupant  are  those  forces  exerted 
by  the  cushions  and  floor,  illustrated  in  figure  11.  The  forces 
of  the  seat  back  pass  through  the  mass  centers  of  the  upper  torso 
and  head  segments,  and  the  normal  forces  applied  to  the  legs  by 
the  seat  bottom  cushion  pass  through  the  mass  centers  of  the  thigh 
segments.  The  contact  surface  for  the  lower  torso  is  not  located 
at  the  mass  center  of  that  segment  but  is  an  ellipsoid  whose  major 
axis  passes  through  both  hip  joints.  The  force  of  the  back  cushion 
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Figure  11.  External  forces  of  cushions  and  floor. 


and  the  normal  component  of  the  seat  bottom  cushion  force  both 
pass  through  the  center  of  this  ellipsoid,  i.e.,  through  a  point 
midway  between  the  hips.  The  normal  components  of  the  floor  and 
footrest  forces  are  applied  to  the  body  at  the  lower  ends  of  the 
leg  segments. 

Each  of  the  forces  mentioned  above  acts  normal  to  the  surface  apply 
ing  the  force,  or,  in  other  words,  its  direction  is  determined  by 
the  plane  of  the  surface,  whether  cushion  or  floor.  As  shown  in 
figure  11,  friction  forces  are  also  applied  by  the  seat  bottom  cush 
ion  and  the  floor.  The  friction  force  is  computed  as  the  product 
of  the  coefficient  multiplied  by  the  normal  force.  It  is  applied 
in  a  direction  opposite  to  the  tangential  component  of  relative 
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velocity  between  the  occupant  segment  and  the  appropriate  cushion 
or  floor  surface.  In  order  to  avoid  abrupt  changes  in  direction 
of  the  friction  force,  the  force  is  reduced  sinusoidally  when  the 
magnitude  of  the  tangential  velocity  drops  below  a  predetermined 
limit.  (A  limiting  value  of  1  ft/sec  has  produced  satisfactory 
results  in  SOM-LA  and  is  thus  used  by  the  program.) 

All  contact  forces  are  calculated  by  first  determining  the  pene¬ 
tration  of  a  contact  surface  on  the  occupant  into  a  surface  with 
known  force-deflection  characteristics.  Using  the  seat  cushion 
force  as  an  example,  the  pertinent  dimensions  of  the  seat  and  the 
parameters  required  to  determine  the  penetration  of  the  abdomino- 
pelvic  segment  (segment  1)  into  the  cushion  are  illustrated  in 
figure  12.  Xp  and  Zp  are  coordinates  of  the  center  of  the  contact 
surface  of  segment  1,  and  R^  is  the  radius  of  the  contact  surface 
in  the  (x^  -  z^)  plane.  (Although  this  contact  surface  is  an  el¬ 
lipsoid,  cross-sections  parallel  to  the  (x^  -  z^)  plane  are  circu¬ 
lar.  The  dimensions  of  the  contact  surfaces  will  be  discussed 
in  section  2.5.)  The  position  of  the  seat  pan  is  defined  by  its 
height  Zg  above  the  origin  of  the  aircraft  coordinate  system  and 
the  angle  ns  that  it  makes  with  the  aircraft  (Xft  -  Yft)  plane. 

The  unloaded  thickness  of  the  seat  cushion  is  t_,  and  the  loaded 
thickness  under  segment  1  is  t.  Summing  the  dimensions  in  the  Zft 
direction  gives 

Zp  «  Zs  +  +  t ) / cos  0g  +  Xp  tan  0S  (37) 

Solving  equation  (37)  for  the  cushion  thickness, 

t  “  (Zp  -  ZQ)  cos  Qq  -  Rj^  -  Xp  sin  0  g  (38) 

The  deflection  of  the  seat  cushion  is  then 

6C  °  te  “  t  (39) 
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Figure  12.  Seat  cushion  deflection. 


and  the  force,  which  is  assumed  to  act  normal  to  the  plane  of  the 
seat  pan  and  pass  through  the  center  of  curvature  of  the  contact 
surface,  is  calculated  from  deflection  according  to 

b6c 

F  =  A(e  -  1)  (40) 

c 

To  each  normal  force,  a  damping  term  is  applied  which  is  propor¬ 
tional  to  the  deflection  rate.  The  damping  coefficient  is  based 
on  the  Rayleigh  formulation  in  which  the  coefficient  is  propor¬ 
tional  to  both  mass  and  stiffness  accordiag  to 

C  =  2  Bm  +  aK  (41) 

For  a  multidegree-of-f reedom  system  there  will  be  a  discrete  damp¬ 
ing  coefficient  associated  with  each  characteristic  mode.  In  a 
continuous  system  there  will  be  an  infinite  number  of  coefficients, 
although  several  modes  will  generally  dominate  the  dynamic  solution. 
It  was  assumed  for  the  formulation  in  SOM-LA  that  the  damping  ratio, 

.  ,  was  constant  for  all  deformation  modes  of  interest.  This  assump¬ 
tion  greatly  simplifies  the  solution  of  equation  (41)  because  the 


32 


damping  coefficient  is  not  dependent  on  the  system  mass.  Equation 
(41)  can  then  be  simplified  to 

C  =  2Ka  (42) 


where  the  stiffness  K  is  the  gradient  of  the  cushion  force- 
deflection  curve,  and  a  is  a  constant  for  the  system. 


The  procedure  used  in  SOM-LA  to  calculate  the  current  value  of  the 
damping  coefficient  is  based  on  the  curent  value  o.c  K,  and  the  con¬ 
stant  a  which  is  determined  from  input  data.  The  slope  of  the  ex¬ 
ponential  load /deflection  curve  at  deflection  5  is 


K 


(43) 


The  user  supplies  a  damping  coefficient,  CQ,  for  the  zero- 
deflection  condition  which  the  program  uses  co  calculate  a  by 
applying  equations  (43)  and  (42),  thereby  resulting  in 

„  Co 

=  2Ko  (44) 


The  constant,  a,  and  the  current  gradient  of  the  cushion  load- 
deflection  curve,  K,  are  used  at  each  time  step  to  determine  the 
instantaneous  cushion  damping  coefficient  from  equation  (42). 

2.4.2  Restraint  System  Forces 

The  method  used  in  calculating  the  forces  exerted  on  the  body  by 
the  restraint  system  differs  considerably  from  that  described  in 
the  preceding  section  for  the  contact  forces.  The  primary  reason 
for  this  difference  is  that  the  restraint  forces  do  not  act  at 
any  fixed  points  on  the  -  ccupant,  but,  rather,  the  points  of  ap¬ 
plication  vary  with  the  restraint  system  geometry. 
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Figure  13.  Restraint  system  configuration  variables. 


Although  other  configurations  can  be  selected  by  the  user,  a  re¬ 
straint  system  consisting  of  a  lap  belt  and  diagonal  shoulder 
strap  will  be  used  as  an  example.  The  restraint  loads  are  trans¬ 
mitted  to  the  occupant  model  through  ellipsoidal  surfaces  fixed  to 
the  upper  and  lower  torso  segments.  These  surfaces  are  shown  in 
figure  13.  The  locations  of  the  anchor  points  A^ ,  A2,  and  A3  are 
determined  by  user  input  along  with  the  webbing  properties.  The 
buckle  B  for  a  single  shoulder  belt  is  located  according  to  an  in¬ 
put  parameter  which  specifies  the  distance  from  the  appropriate 
point,  in  this  case  A^,  along  the  path  of  the  lap  belt.  For  a 
double-strap  shoulder  harness,  the  buckle  is  placed  on  the  abdom¬ 
inal  contact  surface  between  its  intersections  with  the  thigh  sur¬ 
faces. 


for  the  lower  torso,  where 


and 


x22/a22  +  y22/b22  +  z22/c22  *  1  (46) 

for  the  upper  torso,  where 


^2  “  Jj2^2 
c2  *  ^2^2 

and  these  body  dimensions  are  defined  in  section  2.5. 

The  restraint  forces  are  determined  in  the  same  manner  for  both 
the  upper  and  lower  torso.  First,  the  belt  loads  are  calculated 
from  the  displacements  of  the  torso  segments,  and  the  resultant 
force  on  each  segment  is  then  applied  at  the  point  along  the  arc 
of  contact  between  the  belt  and  the  ellipsoidal  surface  where 
the  force  is  normal  to  the  surface. 

Explaining  this  procedure  in  further  detail  for  the  restraint  sys 
tem  configuration  shown  in  figure  13,  for  any  position  of  the  oc¬ 
cupant,  the  coordinates  of  the  left  shoulder,  the  hips,  and  the 
buckle  connection  B  are  calculated  in  the  aircraft  reference 
frame.  The  length  of  each  side  of  the  lap  belt  is  equal  to  the 
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sum  of  the  free  length  in  a  straight  line  from  an  anchor  point 
(A^  or  A2)  to  the  outermost  point  on  the  hip  contact  surface 
(C^  or  C2)  added  to  the  length  of  an  arc  from  that  point  on  the 
hip  to  the  buckle  (B).  The  shoulder  belt  length  is  computed  as 
the  sum  of  the  distance  from  the  anchor  point  ( )  to  a  point  of 
tangency  on  the  top  of  the  shoulder  contact  surface  (D),  the  dis¬ 
tance  from  the  buckle  to  the  extreme  anterior  (forward)  point  on 
the  ellipsoidal  chest  surface  (E),  and  the  length  of  an  arc  over 
the  chest  surface  between  points  D  and  E.  If  the  length  of  the 
belt  segment  should  exceed  the  equilibrium  (zero  load)  length  cal¬ 
culated  initially,  then  there  is  some  tensile  force  in  the  belt. 
The  resultant  force  on  each  segment  is  the  vector  sum  of  the  belt 
forces.  Friction  between  the  shoulder  belt  and  chest  along  the 
length  of  the  belt  is  taken  into  account  by  reducing  the  load  in 
the  belt  between  the  chest  and  buckle  by  a  constant  fraction  of 
the  load  in  the  free  length  between  the  anchor  point  and  the  body 
surface.  The  resultant  force  on  the  lower  or  upper  torso  segment 
may  be  written  generally  as 

P  B  Fx  i  +  Fy  o  +  Fz  k  (47) 

where  Fx,  F^,  and  Fz  are  components  in  the  local,  segment-fixed 
coordinate  system. 

To  find  the  point  on  the  segment  where  F  is  normal  to  the  sur¬ 
face,  consider  first  the  equation  of  an  ellipsoid; 

x2/a2  +  y2/b2  +  z2/c2  =  1  (48) 

which  may  also  be  expressed  in  functional  form  as 

f  (x,  y,  z)  =  x2/a2  +  y2/b2  +  z2/c2  -  1  (49) 

where  the  ellipsoid  can  be  regarded  as  the  level  surface  f=0  of 
the  function.  At  any  point  (x,  y,  z)  on  the  surface,  the  gra¬ 
dient  of  f  is  normal  to  the  surface.  The  gradient  is  given  by 
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(50) 


grad  f  =  (2x/a2)  i  +  (2y/b2)  o  +  (2z/c2)  k 


and  at  the  point  of  application  of  the  resultant  force,  grad  f  is 
col linear  with  F.  Making  use  of  the  proportionality  between  the 
components  of  the  two  vectors, 

Fx  =  Cx/a2 

Fy  =  Cy/b2  (51) 

F  =  Cz/c2 
z 

where  C  is  an  arbitrary  constant.  Solving  equation  (51)  for  the 
coordinates  (x,  y,  z)  and  substituting  into  equation  (48) 


\-C-l  -2  *  [-t-l  -2  [—)  ? 


=  1 


(52) 


which  leads  to 


F  sA2 

X 


/f  b\2  /p  c\2 

+  (-X-j  +  «  1 


C2  =  F  2a2  +  F  2b2  +  F  2c2 
x  y  z 


(53) 


C  =  t  \Fx2a2  ♦  Fy2b2  ♦  Fz2c2 


(54) 


the  point  of  application  of  F  is  then 


x  *  Fxa‘/C 


y  -  FybVC 


z  =  F  c  /C 

Z 
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The  negative  sign  on  C  can  be  explained  by  the  fact  that  each 
coordinate  in  the  local  system  is  opposite  in  sign  to  the  cor¬ 
responding  component  of  the  resultant  force,  or 

x>  0  if  Fx  <  0 

y  >  0  if  Fy  <  0  (56) 

z  >  0  if  F  <  0 
z 

The  capability  of  the  belt's  point  of  application  of  resultant  belt 
loads  to  move  relative  to  the  torso  surfaces  allows  simulation  of 
the  "submarining"  under  the  lap  belt. 

2.5  OCCUPANT  DIMENSIONS  AND  INERTIAL  PROPERTIES 

Characteristics  required  by  the  occupant  model  for  each  of  the 
segments  are  the  length,  mass,  center  of  mass  location,  and  moments 
of  inertia.  Also  required  are  the  axial  and  flexural  stiffnesses 
and  damping  coefficients  for  both  vertebral  elements,  as  well  as 
compliance  characteristics  for  the  chest  and  abdomen.  It  is  as¬ 
sumed  that,  for  each  segment,  a  line  connecting  the  joints  is  a 
principal  axis,  so  that  the  required  moments  of  inertia  are  all 
principal  moments.  For  each  torso  element,  of  the  two-dimensional 
model,  the  center  of  mass  may,  in  general,  be  offset  from  the  spine 
as  shown  in  figure  6.  Moments  of  inertia  are  then  moments  with 
respect  to  axes  located  at  the  mass  center.  Final  data  required 
to  describe  the  occupant  are  radii  of  26  contact  surfaces,  which 
are  ellipsoids  and  spheres. 

For  two  "standard"  occupants,  a  50th-percentile  human  male  and 
a  50th-percentile  anthropomorphic  (Part  572)  dummy,  all  the  rt^ 
quired  data  are  stored  within  the  program.  For  other  nons tarda rt 


,  f 

occupants,  the  above-described? data  must  be  provided  as  input. 

The  human  data  are  based  on  the  U.S.  Air  Force  drawing  board  man¬ 
ikin  (reference  7),  whose  dimensions  were  extracted  by  the  FAA 
Civil  Aeromedical  Institute  and  tabulated  with  appropriate  iner¬ 
tial  properties  in  reference  8.  Dummy  dimensions  and  character¬ 
istics  were  obtained  from  reference  9. 

Most  of  the  dimensions  and  inertial  properties  used  in  SOM- LA  wcro 
taken  from  references  8  and  9,  but  because  fewer  segments  are  used 
in  the  mathematical  model,  some  properties  needed  to  be  combined. 
That  is,  in  the  model,  the  hands  are  combined  with  the  lower  arms, 
the  feet  with  the  lower  legs,  and  the  mid  torso  with  the  lower 
torso.  The  distance  from  the  wrist  pivot  to  the  mass  center  of 
the  hand  was  added  to  the  lower  arm  length;  a  contact  sphere  cen¬ 
tered  on  the  end  of  this  link  accounts  for  the  remaining  reach  of 
the  hand.  Similarly,  the  distance  from  the  ankle  pivot  to  the 
foot  center  of  mass  was  added  to  the  lower  leg  length,  and  a  con¬ 
tact  sphere  was  sized  to  provide  the  correct  distance  between  the 
knee  pivot  and  the  floor.  Composite  mass  center  locations  and 
moments  of  inertia  were  calculated  for  the  lower  torso,  lower 
arms,  and  lower  legs. 

2.5.1  Body  Segment  Dimensions 

The  basic  dimensions  of  the  occupant  segments  that  are  required 
in  writing  the  equations  of  motion  are  illustrated  in  figure  14. 
The  lengths  of  the  segments  are,  in  most  cases,  effective  "link 
lengths"  between  joint  centers,  rather  than  standard  anthropome¬ 
tric  dimensions  based  on  external  measurements.  These  lengths, 
for  the  standard  50th-percentile  occupants,  are  presented  in 
table  1. 

For  segments  other  than  the  torso  segment,  the  distance  of  the 
mass  center  of  segment  n  from  the  end  nearest  the  body  reference 
point  (m^)  is  nn.  The  distance  between  the  mass  center  and  the 
far  end  is  given  by 
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TABLE  1.  BODY  SEGMENT  LENGTHS 

(IN.) 

Segment 

5 Oth- Percentile 

Air crewmember 
(reference  8) 

Part  572 
Dummy 

(reference  9) 

Lower  Torso,  L^ 

9.44 

10.5 

Upper  Torso,  L2 

13.1 

11.5 

Neck,  N 

5.10 

4.88(2) 

Head,  Lg 

8.50(1) 

8.35(2) 

Upper  Arm,  L^ 

11.6 

11.3 

Lower  Arm,  L^ 

14.8(3) 

13.3(3) 

Upper  Leg,  Lg 

17.1 

16.5 

Lower  Leg,  Lg 

18.4<3> 

18.Q(3) 

Spine,  S 

12.4 (3) 

10.85(3) 

Seated  Height 

37.0(3) 

36.0,3> 

(1)  Scaled  from  manikin  drawing. 

(2)  Scaled  from  Part  572  drawing. 

(3)  Calculated. 

Note  that  the  lengths  of  the  tor30  segments,  and  L2,  are  not 
used  in  the  two-dimensional  model  shown  in  figure  5  although 
they  are  used  to  generate  contact  surface  ellipsoids  for  the 
graph!'  -»  display.  However,  the  center  of  mass  distances,  and 
n2,  and  the  spinal  length,  S,  are  used  and 

L1  +  L2  °  P1  *  P2  +  S  (58) 

Also,  the  seated  height  is  equal  to 
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Rl  +  +  S  +  P2  +  N  +  L^/2 


(59) 


where  is  the  radius  of  the  pelvic  contact  ellipsoid  in  the  mid- 
sagittal  plane,  N  is  the  neck  length,  and  is  the  head  segment 
length  (z-direction,  approximately  menton  to  top  of  head,  in  an¬ 
thropometric  terms). 


2.5.2 


Segment  Weights  and  Center  of  Mass  Locations 


Body  segment  weights  and  axial  locations  of  centers  of  mass  for 
the  two  standard  occupants  are  presented  in  table  2.  The  dimen¬ 
sions  e^  and  e2,  by  which  the  torso  mass  centers  are  offset  from 
the  spine  were  determined  to  be  0.20  and  0.70  in.,  respectively, 

based  on  simulation  of  vertical  (+G  )  tests.  A  value  of  1.17  in. 

z 

for  the  corresponding  dimension  for  the  head,  e^»  is  based  on  Part 
572  duramv  drawings. 


2.5.3  Body  Segment  Moments  of  Inertia 

Body  segment  moments  of  inertia  are  presented  in  table  3.  As  in¬ 
dicated,  moments  of  inertia  with  respect  to  lateral  (y)  axes  were 
taken  from  references  8  and  9  for  the  standard  human  and  dummy  oc¬ 
cupants,  respectively.  The  moments  of  inertia  with  respect  to  the 
segment  x-  and  z-axes  were  determined  using  approximations  to  seg¬ 
ment  geometry.  The  torso  and  head  segments  were  approximated  by 
ellipsoids.  Assigning  appropriate  anthropometric  dimensions  to 

the  ellipsoid  axes,  ratios  I  /!.,  and  I  /I„  were  calculated  for 

x  y  z  y 

unit  mass.  These  ratios,  multiplied  by  the  1^  from  reference  8 
or  9,  yave  values  of  I  and  I_  for  the  torso  and  head  segments. 

The  identical  procedure  was  used  for  the  extremities,  except  that 
those  segments  were  approximated  by  solid  circular  cylinders. 
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TABLE  3.  BODY  SEGMENT  MOMENTS  OF  INERTIA  ( lb-in. -sec2 ) 


Segment 

50th- Percentile 
Aircrewmember 
(reference  8) 

Part  572 
Dummy 

(reference  9) 

I 

z 

Lower  Torso 

4.03 

1.19 

0.760 

2.323 

Upper  Torso 

2.37 

3.29 

0.926 

1.70 

Neck 

- 

0.019 

0.0177 

- 

Head 

0.160 

0.199 

0.266 

0.233 

Upper  Arm 

0.131 

0.120 

0.135 

0.022 

Lower  Arm 
(including 

hand) 

0.105 

0.254 

0.185 

0.195 

Upper  Leg 

0.212 

1.41 

1.22 

0.873 

Lower  Leg 
(including 

foot) 

1.28 

1.17 

0.S94 

0.505 

2.5.4  Body  Contact  Surfaces 

Twenty-six  surfaces  are  defined  on  the  body  for  calculation  of 
external  forces  exerted  on  the  occupant  by  the  seat  cushions  or 
restraint  system  and  for  prediction  of  impact  between  the  occu¬ 
pant  and  the  cockpit  interior.  These  surfaces  are  ellipsoids, 
cylinders  and  spheres,  as  shown  in  figure  15.  The  dimensions  of 
these  surfaces,  listed  in  table  4,  were  obtained  from  anthropome¬ 
tric  data  in  references  8  and  9  or  scaled  off  the  drawings  of  the 
manikin  and  dummy.  The  surfaces  and  the  dimensions  required  for 
their  description  are  illustrated  in  detail  in  figure  16. 


Surface  Identification 


1 .  Lower  Torso 

2.  Upper  Torso 

3 .  Head 

4.  Right  Upper  Arm 

5 .  Right  Forearm 

6.  Left  Upper  Arm 

7.  Left  Forearm 

8 .  Right  Thigh 

9.  Right  Lower  Leg 
10-  Left  Thigh 

11.  Left  Lower  Leg 

12.  Right  Knee 

13.  Left  Knee 

14.  Right  Foot 

15.  Left  Foot 

16.  Right  Hip 

17.  Left  Hip 

18.  Right  Shoulder  Joint 

19.  Left  Shoulder  Joint 

20.  Right  Elbow 

21.  Left  Elbow 

22.  Right  Hand 

23.  Left  Hand 

24.  Neck 

25.  Right  Shoulder 

26.  Left  Shoulder 


Figure  15.  Occupant  Contact  Surfaces. 

2.5.5  Joint  Rotation 

The  results  of  several  studies  on  the  limits  of  human  joint  mo¬ 
tion  have  been  published.  Two  of  these  studies,  in  particular, 
were  examined  for  applicability  to  the  occupant  model.  First  of 
all,  Dempster's  (reference  10)  data  on  link  lengths  and  inertial 
properties  were  used,  as  discussed  in  preceding  sections,  so  it  was 
considered  appropriate  to  include  his  joint  data  here.  Glanville 
and  Kreezer  (reference  11)  presented  limits  of  joint  motion  for 
both  voluntary  and  forced  rotation;  their  results  appear,  along 
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TABLE 

4 .  BODY 

CONTACT  SURFACE 

RADII  (IN.) 

Contact 

Surface 

50th-Percentile 
Aircrewmember 
(scaled  from  manikin) 

Part  5  72 

Dummy 

(scaled  from  drawinq) 

Lower  Torso 

4.00 

4.50 

Upper  Torso 

5.00 

4.50 

Neck(1) 

2.00 

2.00 

Head(2) 

3.75 

3.44 

Upper  Arm 

2.10 

1.95 

Lower  Arm 

1.65 

1.85 

Upper  Leg 

3.28 

3.10 

Lower  Leg 

2.23 

2.30 

Hip(3) 

3.56 

3.56 

Shoulder ^  ^ 

2.00 

2.00 

Foot 

1.60 

1.60 

(1)  Neck  circumference  divided  by  2  it. 

(2)  Head  length  (anterior-posterior). 

(3)  Hip  breadth  (sitting)/2  -  hip  link  length  (L„). 

(4)  Shoulder  breadth/2  -  shoulder  link  length  (Lg). 


with  Dempster’s,  in  table  5.  Definitions  of  the  various  joint 
motions  are  illustrated  in  figure  17.  Also  included  in  table  5 
are  the  rotations  required  for  the  Part  572  anthropomorphic  dummy. 

All  of  the  rotations  possible  in  the  mathematical  model  are  in¬ 
cluded  in  table  5  and  figure  17,  but  some  are,  naturally,  more 
important  than  others  in  determining  permissible  ranges  of 
motion  for  the  model.  For  the  head,  ventrif lexion  (8)  is  cer¬ 
tainly  the  most  important  component  of  motion  for  frontal  impact, 
borsiflexion  (A)  may  also  be  important  for  frontal  impact,  but 
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Figure  16.  Body  contact  surface  dimensions 
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•Transverse  plane 


Figure  17.  Motion  diagrams. 

the  angles  reported  are  sufficiently  close  to  those  for  ventr  in¬ 
flexion  to  be  considered  the  same.  Lateral  flexion  (C)  is  cer¬ 
tainly  less  important  since  a  pure  lateral  impact  of  an  aircraft 
would  be  rare  indeed,  and  rotation  (D)  will  have  an  insignficant 
effect  on  model  response.  Therefore,  the  limiting  rotation  (L 


(see  section  2.3)  for  the  neck  joint  (i  =  2)  has  been  taken  as 
the  limit  for  voluntary  ventriflexion,  or  60  degrees.  This  angle 
can  be  found  in  table  6,  along  with  the  limiting  angles  for  the 
other  body  joints.  For  all  of  the  other  angles,  flexion  is  the 
most  significant  component  for  the  type  of  motion  that  can  usually 
be  expected  to  take  place  in  a  crash  environment.  Therefore,  the 
limiting  angles  were  all  taken  as  the  limits  for  voluntary  flex¬ 
ion.  Note  that,  for  the  hip  joint,  the  reference  position  of  the 
body  used  in  the  mathematical  model  includes  90-degree  flexion. 
Therefore,  this  amount  has  been  subtracted  from  the  angle  re¬ 
ported  in  table  5,  which  is  defined  relative  to  the  standard  ana¬ 
tomical  reference  position.  Since  the  seated  position  appears  to 
aid  in  flexion  of  the  hip  joint,  the  largest  angle  in  the  table, 
the  one  given  by  the  Part  572  specifications,  was  used  in  deter¬ 
mining  flc  which  is  thus  given  by  g_  =  120°  -  90°  =  30°. 

/  v 


TABLE 

6.  JOINT  LIMITING  ANGLES 

Joint 

Location 

Angle  - 

i 

(deg) 

1 

Back 

40 

2 

Neck 

60 

3,  5 

Shoulder 

180 

4,  6 

Elbow 

142 

7,  9 

Hip 

30 

8,  10 

Knee 

125 

2.5.6  Body  Stiffnesses 

The  lumbar  spine  and  neck  of  the  two-dimensional  model  possess  ex¬ 
ponential  stiffness  characteristics,  in  the  form  of  equation  (40), 
for  both  axial  and  rational  deformation.  Exponential  stiffness 
characteristics  for  the  abdomen  and  chest  are  used  to  soften  the 
input  force-deflection  values  for  the  lap  belt  and  shoulder  belt, 
respectively. 
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3 . 0  SEAT  MODEL 


The  seat  structure  is  modeled  using  the  finite  element  method  of 
analysis.  This  method  has  been  selected  because  it  is  not  reli¬ 
ant  on  previous  testing,  and  it  has  the  flexibility  to  deal  with 
a  wide  range  of  design  concepts.  The  specific  finite  element  for¬ 
mulation  used  in  the  program  is  based  on  the  WRECKER  II  program, 
developed  at  the  ITT  Research  Institute  (reference  12). 

The  SOM- LA  seat  analysis  includes  triangular  plate  elements, 
three-dimensional  beam  elements,  and  spring  elements.  It  has  the 
capability  to  model  large  displacements,  nonlinear  material  be¬ 
havior,  local  buckling,  and  various  internal  releases  for  beam 
elements.  The  large  displacement  formulation  separates  the  ele¬ 
ment  displacement  field  into  a  rigid  body  rotation  and  transla¬ 
tion  associated  with  a  local  coordinate  system  that  moves  with 
the  element  and  small  element  distortions  relative  to  the  current 
position  of  the  element  coordinate  system.  This  formulation  can 
accommodate  extremely  large  rotations  and  deflections  with  accur¬ 
acy  depending  on  the  size  of  the  elements  relative  to  the  curva¬ 
ture  of  the  structure.  Nonlinear  material  formulation  is  based 
on  a  uniaxial  elastic-plastic  stress-strain  law  for  beam  and 
spring  elements  and  a  biaxial  elastic-plastic  stress-strain  law 
(Von  Mises  yield  criterion)  for  plate  elements.  Internal  releases 
for  beam  elements  include  shear  (transverse  sliding  joint),  moment 
(transverse  hinge  joint),  thrust  (axial  sliding  joint)  and  torque 
(axial  hinge  joint)  releases.  Also,  a  simple  local  buckling  model 
for  thin-walled  tubes  subjected  to  axial  compressive  and/or  bend¬ 
ing  loads  was  incorporated  into  the  program.  This  model  simulates 
the  reduction  in  bending  rigidity  of  the  tube  as  the  cross  section 
distorts  during  local  buckling. 


3.1  SOLUTION  PROCEDURE 


The  solution  procedure  is  based  on  formulation  of  equations  of 
quasi-static  equilibrium  for  the  finite  element  model  in  an  in 
cremental  form  according  to 


”E  if  F? 


— T  <Ui+l  "  “i>  "  \ 


Fi+1  "  Fi 


i+1 


if  F* 

11  1+1 


-  F?)  >  0 


-  F?  <  0 

1  — 


(60) 


where  KT  =  Tangent  stiffness  matrix 

"•“E  thl 

F^  =  External  forces  applied  at  i  solution  time  step 

ui  *  Displacements  (or  rotations)  at  i***1  solution  time  step. 


_ £ 

The  external  forces,  F  ,  including  restraint  system  loads,  occu¬ 
pant  loads  on  the  seat  pan  and  seat  back,  and  seat  support  reac¬ 
tions  are  treated  as  static  loads  on  the  seat  structure.  The 
mass  of  the  seat  structure  is  neglected,  since  in  most  simula¬ 
tions  it  will  be  a  small  percentage  of  the  total  occupant  weight, 

and  the  seat  structure  is  assumed  to  be  in  a  quasi-static  equili- 

— E 

brium  with  the  applied  external  forces  F  . 


During  the  part  of  the  simulation  when  the  external  forces  on  the 
seat  structure  are  increasing,  (  F^^  -  FF  )  >  0,  the  incremental 
displacements,  (Au  B  u^  -  ui\),  are  obtained  from  equation  (60). 
The  total  displacement  of  the  seat  structure  is  obtained  by  sum¬ 
mation  of  the  incremental  displacements.  However,  when  the  ex¬ 
ternal  load(s)  are  decreasing,  (  FF^  -  <_  0,  the  elastic  un¬ 

loading  of  the  seat  structure  is  neglected  and  the  incremental 
displacement  is  assumed  to  be  zero,  (Au  *  0),  since  in  most  sim¬ 
ulations  elastic  deformations  will  be  small  compared  with  plastic 
deformations. 


The  tangential  stiffness  matrix  depends  on  the  state  of  stress 
of  the  seat  structure  and  varies  with  time  during  the  simulation. 
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Therefore,  it  must  be  recomputed,  assembled  and  inverted  at  se¬ 
lected  time  steps  during  the  solution.  The  major  computational 

2 

effort  is  the  inversion  of  K^,  which  requires  NB  multiplications 
where  B  is  the  semibandwidth  and  N  is  the  number  of  degrees  of 
freedom. 

3.2  COORDINATE  SYSTEMS 


Three  different  coordinate  systems  are  used  to  describe  the  finite 
element  model  of  the  seat  structure.  They  are  global,  nodal,  and 
element  coordinate  systems. 

The  global  coordinate  system  (x,  y,  z)  is  fixed  in  space  and 
serves  as  an  inertial  frame  of  reference. 

A  nodal  coordinate  system  (x,  y,  z)  is  attached  to  each  node. 

The  orientation  of  nodal  axes  (x,  y,  z)  with  respect  to  the  global 
axes  at  any  time  is  established  by  the  components  of  three  vectors 
h^,  n2»  n3  which  remain  fixed  along  the  nodal  axes  (x,  y,  z) ,  re¬ 
spectively,  as  the  node  translates  and  rotates.  If  these  three 
unit  vectors  n^,  n2,  h3  form  the  columns  of  a  3  x  3  matrix  N,  then 
any  vector  v  can  be  transformed  from  nodal  to  global  coordinate 
system  by  the  following  time-dependent  transformation 

VG  *  N  vN  {61) 

<■»  #•* 

An  element  coordinate  system  (x,  y,  z)  is  attached  to  each  ele¬ 
ment  and  serves  to  define  the  rigid  body  rotation  and  translation 

A  «  » 

of  the  element.  The  orientation  of  element  axes  (x,  y,  z)  with 
respect  to  the  global  axcrt  at  any  time  is  established  by  the  com¬ 
ponents  of  three  vectors  e^,  e‘2,  which  remain  fixed  along  the 
element  axes  (x,  y,  z),  respectively.,  as  the  element  translates 
and  rotates.  If  these  three  unit  vectors  6^,  e‘2,  Q-j  form  the 
columns  of  a  3  x  3  matrix  E,  then  any  vector  v  can  be  transformed 
from  element  to  global  coordinate  system  by  the  following  time 
dependent  transformation: 
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VG  *  E  Vg  (62) 

A 

For  triangular  plate  elements,  the  z-axis  is  defined  by  the  normal 

r\ 

to  the  plane  formed  by  the  three  corner  nodes,  and  the  x-axis,  by 
a  line  bisecting  the  angle  at  a  selected  node  (1),  as  shown  in 
figure  18. 


.-A 


Deformed 


Undeformed 


Figure  18.  Three-dimensional  plate  element 
(from  reference  12). 


For  beam  elements  (figure  19)  the  x-axis  is  defined  by  a  line 
connecting  the  end  points  of  the  beam,  and  the  y-axis,  by  a  line 

A. 

normal  to  the  x-axis  and  lying  in  a  plane  containing  both  the 
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Deformed 


Unde formed 


Figure  19.  Three-dimensional  beam  element 
(from  reference  12). 


coordinate  reference  point  and  x-axis.  The  remaining  z-axis  is 

A  A 

determined  as  a  normal  to  x-  and  y-axes  by  the  right-hand  rule. 


For  spring  elements,  only  one  element  axis  is  required,  and  it  is 
defined  by  a  line  joining  the  end  points  of  the  spring. 


3.3  ELEMENT  FORMULATION 


Large  displacement  formulation  separates  the  element  displace¬ 
ment  field  into  rigid  body  rotation  and  translation  and  small  ele 
ment  distortion  relative  to  the  current  position  of  the  element 


coordinate  system.  After  the  rigid  body  motion  is  removed,  it  is 
possible  to  use  the  classical  small  deformation  finite  element  for¬ 
mulations.  Consequently,  extremely  large  rotations  and  transla¬ 
tions  can  be  accommodated  with  accuracy  depending  on  the  size  of 
the  elements  relative  to  the  deformed  curvature  of  the  structure. 

3.3.1  Beam  Element 


The  beam  element  is  based  on  the  conventional  small -deflection 
formulation  involving  cubic  displacement  fields  for  transverse 
displacements  and  linear  displacements  for  axial  and  torsional 
displacements. 

From  the  principle  of  virtual  work,  the  general  form  of  the  beam 
element  tangent  stiffness  matrix  (reference  12)  is  given  by 

K-  =  /  DT  C  D  dV  (63) 

1  V 

where  j<T  =>  Element  tangent  stiffness  matrix 

D  =  Matrix  that  relates  element  strains  to  nodal  dis¬ 
placements 

C  *»  Constitutive  matrix  that  relates  stresses  and  strains. 

For  linear  elastic  beams,  the  form  of  the  equation  given  in  equa¬ 
tion  (63)  results  in  the  classical  12  x  12  elastic  beam  stiffness 
matrix  shown  in  figure  20. 

where  A  ®  area  of  cross  section 

E  =  modulus  of  elasticity 
1  ®  moment  of  inertia 
G  »  modulus  of  rigidity 
J  *  torsion  constant 
L  a  length  of  member 
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=  moment  about  x,  (y),  (z)  at  node  i 


force  in  x,  (y),  (z)  direction  at 
node  i 

rotation  about  x,  (y),  (z)  axis  at 
node  i 

x,  (y),  (z)  component  of  displacement 
at  node  i 


For  nonlinear  materials  the  constitutive  matrix  C  is  given  by 

r  -  M 
-  Ae 


(64) 


where  Act 

Af 


incremental  stress 
incremental  strain. 


Using  the  cubic  shape  functions  for  the  transverse  displacements 
and  linear  shape  functions  for  axial  and  torsional  displacements 
in  matrix  D  it  can  be  shown  that  (reference  12) 


K  * 
i±r{' 


Ct?  ' 


'V£2  6yt£j 

Vfx‘ 

Vl  -Sxf2 

‘.*3 

-4xf3 

V? 

-Sx£3 

6y£3 

S  £? 
z  1 

-*x£l 

vi. 

(65) 


where  >  a  element  length 


fj,  »  6x  -  2  St 


f2  ■  6x  -  4* 
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••  v.  V 


f 3  =  f lf 2 
<x  =  /  dA 


A  a 


Ao 


"v  “  '  y  u  dA'  8z  ’  '  *  fi  " 


-  /  yz  if  dA,  s  =  /  y2  if  dA, 


s_  =  /  z2  fa  dA 

2  A  Af- 


(66) 


The  integrals  defined  in  equations  (65)  and  (66)  are  calculated 
numerically  through  the  cross  section  and  along  the  length  of  the 
beam  element.  The  integration  is  piecewise  linear  through  the 
depth  and  linear  along  the  length. 


3.3.2  Spring  Element 


The  spring  element  is  a  one-dimensione-l  element,  and  its  tangent 
stiffness  matrix  is  a  degenerate  case  of  that  of  the  beam  element 
described  in  section  3.3.1. 


The  element  stiffness  matrix  is  given  by 


(67) 


where  k  is  a  spring  constant  that  depends  on  the  state  of  stress. 


3.3.3  Plate  Element 

Triangular  plate  element  formulation  is  based  on  small-deflection 
linear  plate  theory  involving  a  linear  displacement  field  for  mid¬ 
plane  deformations  and  a  cubic  displacement  field  for  plate  bend¬ 
ing  deformations  (reference  13). 
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Node  point  degrees  of  freedom  are  two  in-plane  displacements  u 


and  u  ,  and  two  rotations  6  and  0 

y  x  y 

the  plane  of  the  plate,  0  , 

z 


Nodal  rotations  normal  to 
are  not  admitted.  The  out-of-plane 


displacements  u  are  zero,  since  the  x-y  element  coordinate  plane 
z 

is  established  by  the  deformed  position  of  the  nodes. 

The  development  of  the  tangent  stiffness  matrix  for  plates  is 
similar  to  the  development  for  beams.  The  general  form  of  the 
plate  element  stiffness  matrix  is  the  same  as  equation  (63) 


*T  * 


y~T  -  - 


dV 


where  C  =  constitutive  matrix  that  relates  the  biaxial  stresses 
and  strains 


Aoxx  Ao^  Aqxy 


Ac xx  Aeyy  Aexy 


D  =  Matrix  that  relates  element  strains  to  nodal  dis¬ 
placements  . 

Triangular  plate  element  elastic  stiffness  matrix  is  given  in 
reference  13. 

For  nonlinear  materials,  the  stiffness  matrix  must  be  calculated 
by  numerical  integration  based  on  the  current  state  of  stress  and 
the  plate  strain-displacement  relations.  This  development  given 
in  reference  12  is  as  follows: 


The  plate  strains  can  be  written  as  a  combination  of  the  membrane 
strains  r®  ■  and  curvature  functions  k.  .  in  the  form 


13 


13 


(68) 
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where  u  =  in-plane  displacements 
O’  =  plate  corner  rotations 

j>  =  matrix  of  membrane  strain  functions  (reference  13) 
j;  =  matrix  of  curvature  functions  (reference  12). 


The  complete  strain-displacement  transformation  matrix  D  has  the 
form 


D 


&  0 

0  -z^ 


(69) 


when  this  form  of  D  is  introduced  into  the  stiffness  matrix  in¬ 
tegral  equation  (63)/  it  results  in 


where 


K  m  ^uu  — uO 

— Ou  -00 
•  » 

\u  -  /  iT  c  i  dV 
£ua  •  /  s  /  c  ^  av 

«0U  ■  '  *  /  £  i  dv 
i00  •  r  *2  iT  £  i  «v 


(70) 


(71) 


These  may  be  further  written  in  terms  of  area  and  thickness  in¬ 
tegrals  as 
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K  ■  /  d>  C  <bT  dA 
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(72) 


(73) 


The  matrices  Cuu,  CuQ,  and  CQ6  in  equation  (73)  are  integrated 
numerically  by  trapezoidal  rule.  For  elastic  materials  van¬ 
ishes  and  hence 


The  integrals  in  equation  (72)  are  then  evaluated  by  numerical 
quadratures.  Those  integrations  involve  cross  section  integra¬ 
tions  carried  out  each  time  the  stiffness  matrix  is  created. 


3.4  INTERNAL  RELEASES  IN  BEAM  ELEMENTS 


Internal  releases  for  beam  elements  include  shear  (transverse  slid 
ing  joint),  moment  (transverse  hinge  joint),  thrust  (axial  sliding 
joint)  and  torque  (axial  hinge  joint)  releases.  Combinations  of 
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these  releases  are  possible,  e.g.,  a  ball  joint  is  formed  by  two 
transverse  hinge  joints  in  two  different  directions.  Depending 
on  the  type  of  internal  release,  the  original  beam  element  stiff¬ 
ness  matrix  is  modified  for  force  and  moment  releases  at  speci¬ 
fied  locations  within  the  element. 

For  axial  sliding  and  hinge  joints,  the  original  joint  stiffness 
is  ignored.  For  other  types  of  releases  stiffness  matrix  modifi¬ 
cations,  formulated  in  reference  12,  are  developed  below.  Refer¬ 
ring  to  figure  21,  let  the  release  be  located  at  a  distance  "a" 
from  the  ifc^  end,  and  "b"  from  the  j^*1  end.  It  is  assumed  that 
the  original  stiffness  matrix  K  relating  end  rotations  (<{k,  <j>j) 
and  moments  (M^,  )  of  the  element  is  known.  The  moment  Mg  and 

shear  VD  at  the  joint  are  shown  in  figure  21. 


Figure  21.  Beam  element  internal  releases 
(from  reference  12). 


To  satisfy  the  equilibrium  on  parts  of  the  beam  to  the  left  and 
right  of  the  joint,  a  relationship  between  end  moments  *  (Hi 
and  the  forces  $0  =  (Mg  Vg)  at  the  joint  is  developed 
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(74) 


Sij  =  -  SD 


where 


-1  (-l)na 
1  (-l)nb 


and 


n  = 


1  if  S^j  is  in  z-direction 

2  if  is  in  y-direction 


The  relationship  between  the  forces  and  the  displacements 
Zq  =  (0D,  vD)  will  be 


SD  “  lb  61 


(75) 


Let  KD  be  of  the  form 


— D  = 


D,  ■  D.  . 

ii  ID 


°3i 


33 


It  can  bo  shown  that 


Kd  *  T’1  K  (T*1)7 


(76) 


The  rovers©  transfo motion  also  holds  that 

K  *  T  Kp  TT 


(77) 


At  a  transverse  hinge  )oint  no  moment  can  be  transmitted  about 
the  released  axis*  hence  Mp  *  0.  Substituting  Hp  »  0  into  equa¬ 
tion  (75),  the  reduced  stiffness  matrix  Kp  is 


— D 


D. .  -  D. .  D, . /D. . 

3  j  JO  )i  U 


(78) 
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Then  the  modified  beam  element  stiffness  matrix  can  be  obtained 
from  equation  (77) 

K*  =  T  K£  TT  (78a) 

At  a  transverse  sliding  joint  no  shear  can  be  transmitted  along  ti 
released  axis,  hence  VD  =  0.  Substituting  Vy  =  0  into  equation  (' 
the  reduced  stiffness  matrix  K'  is 


*D  = 


D.  .  -  D. .  D../D.  .  0 

13  31  33 


(79) 


Equation  (73a)  is  also  valid  for  this  case. 


In  addition  to  the  stiffness  matrix,  beam  element  end  moments 
must  also  be  modified.  The  shear  forces  and  moments  at  the  joint 
can  be  expressed  in  terms  of  end  shear  forces  and  moments  by 
solving  for  Sy  from  equation  (74). 


«*  t"1s 


ij 


(80) 


To  implement  a  joint,  apply  -My  to  the  member  at  the  joint,  set 
VD  «  0  but  allow  ftp.  Solvo  equation  (75)  for  0y  and  vQ  in  terras 
of  My.  Solvo  equation  (75)  for  My  in  terras  of  M^  and  .  Com¬ 
bining  the  results  with  equation  (74  ,  the  effect  of  the  release 
on  Sjj  is  given  by 
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To  implement  a  sliding  joint,  apply  -Vy  to  the  member  at  the 
joint,  set  rjy  «  0  but  allow  Vy  and  use  the  same  procedure  as  for 
a  hinge.  Modified  end  forces  take  the  form 
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3.5  MATERIAL  NONLINEARITIES 

Nonlinear  material  formulation  is  based  on  uniaxial  elastic- 
plastic  stress-strain  law  for  beam  and  spring  elements  and  bi¬ 
axial  elastic-plastic  stress-strain  law  (Von  Mises  yield  cri¬ 
terion)  for  plate  elements. 

The  computational  procedure  for  *  biaxial  stress  field  for  plates 
is  presented  in  reference  12  and  is  based  on  the  work  of  Hartzman 
and  Hutchinson,  (reference  14),  as  specialized  for  small  strain  and 
plane  stress  conditions.  The  current  stress  state  at  a  point  in 
the  plate  is  established  as  follows?  Let  the  prior  stress  state 
of  the  point  under  consideration  be  a  3*  0  0  0  and  a  small  in- 

croment  in  strain  from  the  prior  state  to  the  current  state  be 
A.  .  »  At-  First  a  tentative,  current  stress  state,  ci  , 

aa  yy  xy  XX 

11  „ „  is  calculated  as  though  the  strain  increment  were  com- 

yy 

pletely  elastic. 


0  0 
0  O  a  w  + 

XX  XX 


~  (Ae  ♦  vAc  ) 

,2  yy  yy' 


e  0 

O  o  0  v  ♦ 

yy  yy 


(Ar.  ♦  vAc  ) 

yy  xx' 


«i  ♦  G  Ac 

xy  xy 


where  L  -  Modulus  of  Elasticity 
G  0  Shear  Modulus 


v  =  Poisson’s  Ratio 
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An  effective  stress,  o&,  to  determine  whether  plastic  flow  has 
taken  place  during  the  s crain  increment  is  calculated  using  the 
Von  Mises  criterion 


(5, 


XX 


35  2>1/2 
xy ' 


(84) 


If  o  is  less  than  the  prior  effective  stress  0e^  at  which  yield¬ 
ing  occurred,  elemental  loads  are  decreasing  and  the  tentative 

stresses  calculated  (o  ,  d ,  a  }  are  the  correct  values  of  the 

xx  yy  xy  ^ 

current  stress  values  (c-  ,  a.  a).  If  o_  is  greater  than  the 

xx  yy  xy  e 

prior  effective  stress  a®,  at  which  yielding  occurred,  tentative 
stresses  calculated  must  be  modified  to  account  for  the  plastic 
behavior.  Hartsman  and  Hutchinson  (reference  14)  have  shown  that 
the  true  value  of  the  current  effective  stress  will  be 
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where 


(E)(E  ,) 

Hi  “ 

pi 


i  a  1,  2 


(89) 


and 


F  ^  «  Plastic  modulus  (figure  22) 


The  current  state  of  stre-s  is  then  given  by 
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where  *  is  given  by 


1  [  °e 

5  K 
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(87) 
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3.6  LOCAL  BUCKLING  OF  BEAM  ELEMENTS 


Local  buckling  is  one  of  the  failure  modes  of  thin-walled  tubes, 
typically  used  in  light  aircraft  seats,  when  subjected  to  axial 
compressive  and/or  bending  loads.  Evaluation  of  the  results  of 
dynamic  tests  of  seat-occupant  systems  have  indicated  that  the 
local  buckling  of  thin-walled  tubes  have  a  significant  effect  on 
the  response  of  the  seat-occupant  system.  Therefore,  a  simple 
local  buckling  model  for  thin-walled  tubes  was  incorporated  into 
the  program. 

Local  buckling  is  traditionally  expressed  in  terms  of  a  moment¬ 
bending  curvature  diagram  as  shown  in  figure  23. 


Curvature  change  -  6  Buckled 

Figure  23.  Moment  capability  versus  curvature 
of  a  thin-walled  circular  tube. 


od 


The  cross  section  goes  through  several  stages  of  deformation  as 
the  structure  bends.  Although  during  this  cross-section  distor¬ 
tion  axial  stresses  redistribute  themselves,  it  was  proposed  in 
reference  15  that  the  reduced  bending  rigidity  is  most  strongly 
related  to  the  loss  of  lateral  moment  arm  of  the  axial  forces. 
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For  a  circular  tube  under  uniform  compression  the  stress,  cr  ,  at 

Li 

which  local  buckling  starts  is  given  in  reference  16  as 

°L  *  kce(5>  <88> 

where  0.4  <  K  <1.2  (recommended  values) 
c 

E  =  modulus  of  elasticity 
t  =  wall  thickness  of  the  tube 
D  =  diameter  of  the  tube. 

It  was  also  proposed  in  reference  16  that  for  tubes  subjected  to 
bending,  the  local  buckling  stress,  o  ,  as  a  function  of  material 
and  geometry,  can  be  taken  as  that  for  a  circular  tube  under  uni¬ 
form  axial  compression,  since  in  bending,  a  significant  portion 
of  the  circumference  is  subjected  to  a  relatively  uniform  com¬ 
pression  field.  However,  dynamic  tests  of  seat-occupant  systems 
have  indicated  that  the  thin-walled  circular  tubes  have  suffered 
local  buckling  at  stress  levels  much  below  that  predicted  by 
equation  (88).  This  was  partly  due  to  the  fact  that  the  bending 
stresses  have  exceeded  the  yield  point  and  the  tubes  have  also 
suffered  plastic  deformation. 

To  account  for  plasticity  as  well  as  the  cross-sectional  proper¬ 
ties,  an  empirical  relation,  equation  (89)  was  used  in  the  pro¬ 
gram  to  predict  the  local-buckling  stress 

=  °y  '  <89) 

a  local  buckling  stress 
=  yield  stress 

»  local  buckling  coefficient 


where 


y 

K 
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t  =  wall  thickness  of  the  tube 
D  =  diameter  of  the  tube*. 

The  cross  sections  of  beam  elements  are  defined  as  thin-walled 
plate  segments  by  specifying  their  end  points  and  thickness  as 
part  of  the  input  data  (figure  24).  The  tangent  stiffness  matrix 
for  nonlinear  beam  elements  is  then  computed  by  evaluating  the 
integrals  given  in  equation  (66)  by  Gaussian  quadrature  over  the 
plate  segments  defining  the  beam  cross  section. 


Z 


Figure  24.  Circular  tube  cross  section  defined 
by  eight  thin-walled  plate  segments. 

The  program  computes  the  stresses  at  all  Gaussian  stations  across 
the  cross  section  at  each  end  of  the  beam  elements.  These  stresses 
are  then  compared  with  the  local  buckling  stress,  computed 

from  equation  (89)  for  each  beam  element.  If  any  of  the  compres¬ 
sive  stresses  at  the  Gaussian  stations  exceed  a L,  the  deformation 


* por  rectangular  tube s  the  diameter  of  a  circumscribing  circular 
tube  is  used. 
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of  the  cross  section  is  modeled  by  modifying  the  radial  location 
of  that  Gaussian  station  using  the  following  expression 


ai+l I  >  aL 


(90) 


where  R.  =  radial  location  of  the  Gaussian  station  at  time 
1  step  i 

=  local  buckling  stress 

°i  +  l  =  comPress;i-ve  stress  at  the  Gaussian  station  at 
time  step  i+1 

K  *  local  buckling  constant  (0.50  recommended). 


Consequently,  reduced  bending  rigidity  of  the  cross  section  due 
to  the  decrease  in  the  lateral  moment  arm  of  axial  forces  during 
local  buckling  can  be  modeled. 


4 . 0  SIMULATION  COMPUTER  PROGRAM 


The  digital  computer  program  based  on  the  occupant  and  seat  models 
described  in  chapters  2  and  3  is  called  Seat/Occupant  Model  -  Light 
Aircraft  (SOM-LA) .  It  has  been  written  entirely  in  FORTRAN  to 
ensure  a  high  degree  of  compatibility  with  various  digital  com¬ 
puter  systems.  During  development,  the  program  has  been  run  on 
IBM,  UNIVAC,  and  CDC  computer  systems. 

The  elements  of  the  program  can  be  considered  in  terms  of  three 
general  operations: 

•  Input  and  initialization 

•  Solution 

•  Output 

which  are  summarized  below  and  discussed  in  detail  in  the  sec¬ 
tions  following.  The  general  flow  of  the  program  is  illustrated 
in  figure  25.  Input  data  describing  the  occupant  and  crash  condi¬ 
tions  are  read  first.  If  the  user  requests  output  of  the  predic¬ 
tion  of  impact  between  the  occupant  and  the  aircraft  interior, 
the  coordinates  defining  the  cockpit  surfaces  are  read.  Finally, 
the  seat  data,  either  simple  dimensions  describing  a  rigid  seat 
model  or  detailed  design  data  for  the  finite  element  seat  analysis, 
are  provided.  Based  on  the  input  data,  the  values  of  constants, 
such  as  occupant  dimensions  and  properties  are  calculated,  and  the 
initial  position  of  the  occupant  is  determined. 

The  solution  loop  is  entered  for  the  first  time  with  the  aircraft 
initial  velocity  and  the  occupant  initial  position.  At  each  sub¬ 
sequent  entrance  to  the  loop,  the  current  aircraft  displacement, 
velocity,  and  acceleration  components  are  calculated.  The  equa¬ 
tions  of  motion  for  the  occupant  are  set  up  and  solved.  If  a 
finite  element  seat  model  is  being  used,  the  forces  applied  to 
the  seat,  such  as  the  cushion  forces,  are  provided  to  the  seat 
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Figure  25.  Program  flow  chart 


routines  for  computation  of  seat  displacements.  At  time  incre¬ 
ments  equal  to  a  predetermined  print  interval,  the  output  vari¬ 
ables  requested  by  the  user  are  stored  for  printing  after  comple¬ 
tion  of  the  solution. 

4.1  PROGRAM  INPUT 

Input  data  are  read  by  the  program  in  the  following  seven  blocks: 

1.  Simulation  and  output  control  information 

2.  Cockpit  description  (optional) 

3.  Cushion  properties 

4.  Restraint  system  description 

5.  Crash  conditions 

6.  Occupant  description 

7.  Seat  design  information. 

4.1.1  Simulation  Control  Information 


The  first  block  of  data  contains  the  information  required  for 
controlling  execution  of  the  program.  The  initial  time  step  for 
integration  of  the  equations  of  motion,  the  total  length  of  sim¬ 
ulation,  the  number  of  cases  to  be  run,  the  system  of  units  (SI 
or  English),  selection  of  two-or  three-dimensional  occupant,  and 
identification  of  the  desired  output  are  provided  here. 

4.1.2  Cockpit  Description 

For  prediction  of  impact  between  the  occupant  and  the  cockpit 
interior,  ten  plane  surfaces  are  used  to  describe  the  cockpit. 

As  shown  in  figure  26,  six  of  these  surfaces  are  normal  to  the 
Xft  -  ZA  plane  and  four  are  normal  to  the  plane.  The 

first  five  planes  can  be  used  to  describe  the  environment  of  a 
crewseat,  in  which  case  they  represent  the  firewall,  instrument 
panel,  and  windscreen,  or,  for  analysis  of  a  passenger  seat, 
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environment 


Crewseat  environment 
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Figure  26.  Planar  surface  approximation 
to  aircraft  interior. 


they  can  be  rearranged  to  describe  a  seat  back.  Input  data  in¬ 
clude  X  and  2  coordinates  to  define  planes  1-6  and  Y  and  2  coor¬ 
dinates  for  planes  7-10. 

4.1.3  Cushion  Description 

Cushion  load-deflection  characteristics  are  described  by  an  ex¬ 
ponential  function,  whose  coefficients  are  provided  as  input  data. 
The  equilibrium  (zero  load)  thickness  for  bovh  the  seat  and  back 
cushions  are  also  given.  The  cushion  damping  coefficient  for 
zero  deflection  described  in  section  2.4.1  is  also  entered. 


4.1,4  Restraint  System  Description 


The  restraint  system  used  in  the  simulation  may  consist  of  a  lap 
belt  alone  or  combined  with  a  single-  or  double-strap  shoulder 
harness.  A  lap  belt  tiedown  (negative  -G)  strap  can  also  be  in¬ 
cluded.  The  webbing  force-elongation  curve  is  approximated  by 
three  linear  segments,  which  are  described  by  input  of  points  on 
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the  curve.  The  force  is  computed  by  linear  interpolation  in  this 
table,  as  described  in  section  4.3.1.  The  slack  in  the  webbing 
is  also  provided  by  input  in  units  of  length. 

The  anchor  points  for  the  lap  belt,  shoulder  harness,  and  tie¬ 
down  strap  are  located  by  input  of  rectangular  coordinates  in  the 
aircraft  reference  system.  For  a  double-strap  shoulder  harness, 
the  buckle,  or  point  of  connection  to  the  lap  belt,  is  assumed 
located  on  the  mid-point  of  the  lap  belt.  If  a  tiedown  strap 
is  included  in  this  case  it  also  connects  to  the  other  restraint 
system  components  at  the  buckle.  For  a  single  shoulder  belt, 
which  may  pass  over  either  the  left  or  the  right  shoulder,  an 
input  parameter  locates  the  buckle  by  the  length  of  webbing  be¬ 
tween  the  buckle  and  the  lap  belt  anchor  point.  This  length  may 
be  zero  if  the  buckle  attaches  directly  to  a  rigid  anchor  point. 

4.1.5  Crash  Conditions 


The  aircraft  crash  conditions  are  defined  by  the  initial  velocity 

and  attitude  and  the  acceleration  as  a  function  of  time.  Six 

components  of  velocity  are  required:  three  translational  in  the 

aircraft  coordinate  system  (Vx  »  vy  '  vz  *  and  ^aw'  P*tc^' 

AAA 

and  roll  rates  {$A,  $A).  Each  of  the  six  acceleration  com¬ 

ponents,  which  define  the  acceleration  of  the  aircraft  coordinate 
system,  is  described  by  sixteen  points  in  time  and  acceleration. 

An  example  of  an  approximation  to  an  actual  acceleration  pulse  is 
illustrated  in  figure  27.  Although  many  of  the  higher  frequency 
oscillations  observed  in  the  actual  pulse  probably  contribute 
little  to  the  overall  response  of  the  occupant,  the  use  of  a  large 
number  of  points  reduces  the  effect  of  the  investigator's  subjec¬ 
tivity  in  the  approximation. 

4.1.6  Occupant  Description 

Because  it. has  been  assumed  that  the  principal  user  of  this  pro¬ 
gram  is  interested  primarily  in  the  seat  or  restraint  system,  a 
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Figure  27.  Piecewise  approximation  to  aircraft 
acceleration  component. 


minimum  of  information  ie  required  to  describe  the  occupant. 

Data  for  standard  human  and  dummy  occupants,  as  described  in  sec¬ 
tion  2.5,  are  stored  within  the  program.  Additional  data  must  be 
provided  for  nonstandard  occupants. 


4.1.7  Seat  Design  Information 

The  input  data  required  to  describe  the  seat  consist  of  nodal  co¬ 
ordinates,  material  properties,  cross-section  geometries,  element 
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locations,  and  attachment  conditions,  as  described  earlier.  For 
possible  use  in  restraint  system  or  cabin  configuration  analyses 
where  detailed  seat  response  may  not  be  important  or  seat  design 
unknown,  a  rigid  seat  model  can  be  selected.  Input  data  for  the 
rigid  seat  option  consist  only  of  locations  of  the  seat  pan  and 
seat  back. 

4.2  OCCUPANT  INITIAL  POSITION 

The  initial  position  of  the  aircraft  occupant  is  computed  from 
the  input  parameters  shown  in  figure  28.  It  is  assumed  that  the 
occupant  is  seated  symmetrically  with  respect  to  the  aircraft 
(X^  -  2ft)  plane  or,  equivalently,  that  the  segment-fixed  yn~axes 
are  all  parallel  to  the  Y^-axis.  The  angular  coordinates 
Y^  (i  =  1,  2,  3,  4)  define  the  rotation  of  segments  1-4  relative 
to  the  axis  and,  because  of  the  symmetry  condition,  segment  6 
is  parallel  to  segment  4.  Positive  angles  are  shown  in  figure  28 
although  Vj  and  Y2  are  usually  negative,  considering  the  torso  to 
be  approximately  parallel  to  an  aftward-sloping  seat  back.  The 
angle  describes  the  position  of  the  forearms  relative  to  the 
upper  arms,  and  is  the  initial  value  of  a,.  and  The  distance 

X(1  is  the  initial  X-coordinate  of  the  heels  (the  inferior  ends  of 
segments  9  and  11).  The  procedure  described  below  consists  of 
seating  the  occupant  in  such  a  position  that  static  equilibrium 
is  achieved  among  the  forces  exerted  by  the  seat  cushion,  floor, 
and  either  the  restraint  system  or  the  back  cushion. 

The  first  step  in  determining  the  initial  position  for  the  three- 
dimensional  occupant  involves  calculating  the  Euler  angles  for  the 
torso,  head,  and  arm  segments,  since  this  procedure  does  not  re¬ 
quire  consideration  of  the  forces  due  to  the  cushions  and  floor. 
Because  the  input  parameters  illustrated  in  figure  28  define  the 
position  of  the  occupant  in  the  aircraft  coordinate  system,  the 
orientation  of  the  aircraft  must  be  described  in  the  inertial  sys¬ 
tem.  For  an  aircraft  in  level  flight  with  zero  pitch,  roll,  and 
yaw,  it  is  assumed  that  the  aircraft  coordinate  axes  (Xft,  Yft,  Z A) 
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Figure  28.  Initial  position  input  parameters. 

are  parallel  to  the  fixed  coordinate  axes  (X,  Y,  2)  at  the  ini¬ 
tial  time  (t  s  0).  A  general  orientation  of  the  aircraft  refer¬ 
ence  frame  is  obtained  by  the  same  sequence  of  rotations  defined 
in  soction  2.1.1  for  the  body  segments.  Defining  the  rotations 

i  Yaw 
A 


Pitch 


the  orientation  of  the  aircraft  relative  to  the  inertial  system 
is  described  by  the  coordinate  transformation 


a  <  y 


where  the  elements  of  (A]  are 


A..  a  cos  <{/.  cos  0, 


A,  7  =  cos  i sin  Qa  sin  -  sin  cos  <p. 


A,  a  cos  i sin  o.  cos  *■  sin  _  sin 


A2i  °  sin  cos 


A2?  “  sin  v*A  sin  sin  ♦  cos  ^  cos 


A7;j  “  sin  sin  0^  cos  ^  -  cop  sin 


ft3l  *  -sin  "A 


*32  .  cos  »A  sin  ?A 


A33  °  COS  ^A  008  ^A 


The  rotation  of  body  segment  n  relative  to  the  aircraft,  remem¬ 
bering  that  the  symmetry  condition  requires  that  yn  is  parallel 
to  Y^,  is  described  by 
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=  A21  cos  yn  -  A23  sin  y 


A21  sin  Yn  *  A23  cos  y, 


A31  cos  Yn  '  A33  sin  Yn 


Bn  =  A 
“32  "32 


B33  =  A31  sin  y„  +  A,,  cos  y 


n  "33 


Comparison  of  equation  (94)  with  equation  (1)  points  out  that  the 
transformation  matrices  [Tn]  and  [Bn]  are  equivalent.  Because 
[Tn]  is  a  function  of  Euler  angles  for  segment  n,  equating  the 
elements  of  [Tn],  and  [Bn]  through 


[Tn]  *  [Bn] 


n  *  1,  2,  3,  4,  6/8,  10 


permits  calculation  of  the  initial  values  of  the  generalized  co¬ 
ordinates  from  input  parameters  Yn  and  0^,  and  The  pro¬ 
cedure  as  used  in  Program  SOM-LA  is  outlined  below. 


First  0  is  determined  as  follows: 
n 


*31  *  B31  or  _sin  6n 


which  gives 


en  =  sin'1  (-B^) 


The  cosine  is  then  found  by 


cos  0n  =  cos  [sin“*  ( ~B ^ ^ ) 3 
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so  that  t|j  can  be  determined: 

T1} .  =  B1? .  or  cos  i|*  cos  6  =  B “ 

11  11  n  n  11 

gives 

i|»n  «  cos"1  (bJj/cos  9n)  (99) 

and,  for  determination  of  <J>  , 

T33  =  B33  or  cos  °n  cos  =  B33 

gives 

4>n  «  cos"1  (B33/COS  0R)  (100) 

Equations  (98)  through  (100)  are  used  for  segments  1,  2,  3,  and  4} 
the  symmetry  requirement  provides  the  Euler  anglos  for  segment  6, 
which  are  equal  to  those  for  segment  4.  At  this  point  the  gen¬ 
eralized  coordinates  q4  through  q2Q  have  been  determined.  The 
next  step  involves  seating  the  occupant  and  calculating  Xj ,  Y^, 
and  (q^,  q2,  and  q^)  from  static  equilibrium.* 

Because  the  problem  of  seating  the  occupant  is  statically  inde¬ 
terminate,  certain  simplifying  assumptions  arc  made.  The  first 
assumption,  which  is  approximately  correct  for  typical  seating 
positions,  is  that  15  percent  of  the  occupant's  weight  is  sup¬ 
ported  by  the  floor.  In  other  words,  85  percent  is  supported  by 
the  seat  cushion  and,  depending  on  the  aircraft  attitude,  the 
restraint  system  or  the  back  cushion. 


*Note  that  the  computation  of  Euler  angles  is  required  only  for 
the  three-dimensional  occupant  model.  For  the  two-dimensional 
model,  the  corresponding  generalized  coordinates  are  obtained 
directly  from  input  data. 
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A  first  approximation  to  the  initial  position  is  made  for  the  as¬ 
sumption  of  level  flight  (0A  =  =  <J>A  =  0).  The  cushion  forces 

act  on  the  body  as  shown  in  figure  29,  where  it  is  assumed  that 
15  percent  of  the  occupant  weight  is  supported  by  the  floor,  as 
discussed  in  the  preceding  paragraph.  Summing  forces  gives 


Fx  :  F0  cos  eB  -  Fg  sin  6g  =  0 
A  A 


Fz  :  Fg  sin  0g  +  Fg  cos  6g  =  W’ 


which  can  be  solved  for  the  cushion  forces: 


Fg  =  W'  cos  eB/cos  (eB  -  0g) 


FB  *  W'  sin  0g/cos  (0g  -  0g) 


(101) 


(102) 


Dimensional  considerations  permit  the  coordinates  of  point  P  to 
be  written  as  functions  of  the  thicknesses  tg  and  tg  of  the  com¬ 
pressed  seat  and  back  cushions,  respectively. 


Zp  =  Zg  +  (Rb  +  tg ) / cos  0g  +  (Xp  -  Xg)  tan  0g 


Xp  =  Xs  +  (Rb  +  tg ) / cos  Og  -  (Zp  -  Zg)  tan  0fi 


which  can  be  solved  for  X„  and  Z_  to  give 

P  P 


Xp  =  Xg  -  fB  sin  0 g  +  f2  cos  0g 


Zp  =  Zg  +  f^  cos  0 g  +  f2  sin  0g 


where 


fi  «  (Rj_  +  tg ) / COS  ( 0g  -  0  g) 


(103) 


(104) 


f 2  “  (R1  *  V/cos  teB  "  es) 


85 


Figure  29.  Forces  acting  on  occupant 
torso  (level  flight). 

Since  the  force-deformation  characteristics  of  the  cushions  are 
known  from  input  data,  the  compressed  thicknesses  tg  and  tg  can 
be  calculated  from  equation  (102).  These  values,  when  used  in 
equation  (104),  give  the  coordinates  of  point  P  for  the  first  ap¬ 
proximation  of  level  flight.  The  equilibrium  (zero-load)  lengths 
of  the  lap  belt  and  shoulder  belt(s)  are  calculated  for  the  body 
in  this  position. 

Next,  the  aircraft  is  rotated  to  the  attitude  specified  by  the 
input  conditions  of  pitch,  roll,  and  yaw.  Nose-up  pitch  will 
tend  to  load  the  back  cushion,  and  the  analysis  will  be  the  same 
as  that  described  above  for  level  flight,  except  that  the  W  vec¬ 
tor  in  figure  29  will  have  a  component  in  the  XA  direction. 
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Nose-down  pitch,  on  the  other  hand,  will  tend  to  load  the  re¬ 
straint  system.  An  iterative  procedure  is  used  to  determine  the 
correct  position  for  this  case.  Referring  to  figure  30,  summing 


where  0  is  the  angle  between  the  floor  and  the  projection  of  the 
lap  belt  on  the  Xft  -  ZA  plane.  An  initial  estimate  to  0^  is  made 
using  the  body  position  calculated  by  equation  (104)  for  the  level 
flight  assumption.  The  angle  is  defined  according  to 

0L  =  sin"1  [ (2p  -  ZL)/  ^(Xp  -  XL)2  +  (Zp  -  ZL)2]  (106) 

The  forces  in  the  seat  cushion  F-  and  the  lap  belt  F,  are  deter- 

o  Li 

mined  using  this  value  of  0L  in  equation  (105).  From  the  input 
force-deformation  characteristics  for  the  seat  cushion  and  lap 
belt,  the  deformations  6g  and  <5^  are  calculated.  These  deforma¬ 
tions  are  used  to  determine  new  values  of  Xp  and  Zp?  this  proce¬ 
dure  amounts  to  permitting  the  body  to  further  compress  the  seat 
cushion  and  slide  forward  into  the  lap  belt.  Following  through 
the  procedure,  the  new  length  for  one  side  of  the  lap  belt  is 


LL  =  LLe  +  6L 


(107) 


where  LLe  is  the  equilibrium  length.  The  new  value  of  Xp  is 
given  by 


Xp  *  XL  ♦  (<Ll  -  Lh)2  -  (Yp  -  YL)2]  cos  0l  (108) 

where  is  one-half  the  hip  breadth  and  Yp  is  the  Y-coordinate 
of  the  right  hip  in  the  aircraft  system.  Tho  new  value  of  Zp  is 
computed  for  the  new  cushion  thickness  tg  using  equation  (103), 
which  is  repeated  here  for  continuity: 

Zp  =  Zs  +  (R1  +  tg ) / cos  0S  ♦  (Xp  -  Xg)  tan  0g 

The  new  occupant  position,  determined  by  •jqt’stions  (108)  and  (103) 
is  used  in  equation  (106)  to  recalculate  the  lap  belt  angle 
and  the  procedure  is  repeated  until  two  consecutive  values  of  Xp 
differ  by  less  than  5  percent.  The  coordinates  of  the  mass  center 
of  segment  1  (X^,  Y^,  2 ^ )  are  then  calculated  from  Xp,  Yp,  and  Zp. 
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At  this  point  the  generalized  coordinates  for  the  upper  body  have 
been  determined.  The  final  task  will  be  to  determine  the  coor¬ 
dinates  for  the  legs.  Referring  to  figure  31,  the  angles  y8  and 
can  be  found  from  simple  geometric  relationships  among  the  di¬ 
mensions  shown.  The  Euler  angles  4»g/  6g»  and  $g  are  obtained 
from  y8  by  using  equations  (98)  through  (100),  and  the  correspond¬ 
ing  coordinates  for  segment  10,  by  symmetry.  The  knee  angles  are 
given  by 

am  =  11  “  °k  »  ■  9,  U  (109) 

to  complete  the  initialization  of  the  generalized  coordinates. 


Figure  31.  Leg  position. 

4.3  PROGRAM  SOLUTION  PROCEDURE 

The  first  operation  in  each  solution  step  includes  the  calcula¬ 
tion  of  new  values  for  the  aircraft  acceleration  components  and 
their  subsequent  integration  to  obtain  aircraft  velocity  and 
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displacement  components.  Then  the  matrix  form  of  the  equations 
of  motion,  using  equation  (19) 

(A(q) ]  {q}  =  (B(q,q)}  ♦  {P(q)}  +  (R(q,q))  +  (Q(q,q)} 

are  set  up  for  solution  and  solved,  as  discussed  below. 

4.3.1  Setup  of  Equations  of  notion 

The  elements  of  (A),  {B),  {P},  and  {R}  are  calculated  using  the 
current  values  of  the  generalized  coordinates  and  velocities. 

The  elements  of  (Q),  which  is  the  vector  of  generalized  external 
forces,  are  calculated,  as  discussed  in  section  2.3.  The  exter¬ 
nal  forces  depend  on  displacements  of  the  aircraft,  which  deter¬ 
mine  the  notion  of  the  seat,  floor,  and  restraint  system  anchor 
points  relative  to  the  body.  Prom  these  displacements  new  deflec¬ 
tions  of  the  cushions,  floor,  and  restraint  system  are  calculated. 

4.3.2  Solution  of  Equations  of  Motion 

The  system  of  equations  is  solved  for  the  generalized  accelerations 
by  first  combining  the  vectors  on  the  right-hand  side: 

(A)  { q )  «  (BM  (110) 

where  «BM  ■  (B)  ♦  (P)  ♦  (R)  ♦  {Q) 

and  solving  for  (ql  using  Grout  decomposition  followed  by  Porsythe- 
Molor  elimination. 


The  resulting  set  of  N  second-order  differential  equations  have 
the  general  form 


q3  (t»0)  *  q^0  (t-0)  «  q jQ  j  ■  2,...,  N 


(111) 
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where  N,  the  number  of  degrees  of  freedom,  is  12  for  the  two- 
dimensional  model  and  29  for  the  three-dimensional  model. 


These  equations  can  be  rewritten  as  2N  first-order  equations 
having  the  general  form 


Vj  -  fj  (t,  yx,  y2, ,  yN,  q1#  qN) 

4j  =  Y  j  (U2) 

Vj  (t=0>  -  4j0  qj  (t«0)  =  9 j o 

Numerical  integration  of  this  set  of  equations  is  accomplished, 
using  the  Adams-Moulton  predictor-corrector  method  with  a  vari¬ 
able  step  size.  This  method  uses  the  difference  equations 


»  y.  , ♦  (55f •  -  59f .  „  .  ♦  37f .  „  ,  -9f .  n  „)  (113) 

•'j,n*l  24  j,n  j,n-l  j  ,n-2 

as  the  predictor  and 


yiC)  ay  ♦4“  (9f^  ♦  19f . 

y3,n*l  yj,n  ^4  £j#n*l  j,n 


5£j,n-l  *  fj,n-2'  <ll«> 


as  the  corrector.  Starting  values  are  provided  by  the  classical 
fourth-order  Uunge-Kutta  method.  Input  data  includes  upper  and 
lower  error  bounds  for  the  solution.  Error  bounds  for  each  vari¬ 
able  are  calculated  and  compared  at  each  step  with  the  difference 
between  the  predicted  value  yjp*  and  the  corrected  value  yjc* .  If 
this  difference  exceeds  the  upper  bound  for  any  j,  the  step  size 
is  halved.  If  this  difference  is  less  than  the  minimum  error 
bound  for  all  j  and  for  three  successive  steps,  the  step  size  is 
doubled.  Halving  the  step  si2e  is  accomplished  by  interpolation 
of  past  data,  whereas  doubling  is  effected  by  alternate  selection 
of  past  data.  The  solution  can  be  run  with  a  fixed  step  size  by 
making  the  upper  and  lower  error  bounds  prohibitively  large  and 
small,  respectively,  or  by  using  equal  values  for  the  maximum  and 
minimum  step  si2e  which  are  also  included  among  input  data. 
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4 . 4  PROGRAM  OUTPUT 


Output  data  consist  of  ten  blocks  of  information  that  are  selected 
for  printing  by  user  input.  The  data  include  time  histories  of 
the  variables,  which  are  simply  sorted  during  the  solution  at 
predetermined  print  intervals  as  follows: 

1.  Occupant  segment  positions  in  the  aircraft  coordinate 
system  (X,  Y,  2,  pitch,  and  roll). 

2.  Occupant  segment  velocities  in  the  aircraft  coordinate 
system  (X,  Y,  and  2). 

3.  Occupant  segment  accelerations  in  the  segment-fixed 
coordinate  systems  (x,  y,  z,  and  resultants). 

4.  Restraint  system  loads. 

5.  Cushion  loads. 

6.  Aircraft  displacement,  velocity,  and  acceleration. 

7.  Injury  criteria. 

8.  Details  of  contact  between  the  occupant  and  the  air¬ 
craft  interior. 

9.  Seat  structure  nodal  forces. 

10.  Seat  structure  element  stresses. 

Printer  plots  are  provided  for  occupant  segment  accelerations,  re¬ 
straint  system  loads,  and  cushion  loads.  The  option  of  two  dif¬ 
ferent  filters  is  also  provided  for  the  occupant  segment  acceler¬ 
ations  and  cushion  loads. 

Two  of  the  above  blocks  of  output  data  will  be  discussed  in  fur¬ 
ther  detail. 

4.4.1  Impact  Prediction 

For  prediction  of  impact  between  the  occupant  and  the  cockpit  in¬ 
terior,  26  surfaces  are  defined  on  the  body.  These  surfaces  were 
illustrated  in  figure  15,  and  their  dimensions  discussed  in  sec¬ 
tion  2.5. 
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The  distance  between  each  of  these  occupant  contact  surfaces  and 
the  aircraft  cockpit  surfaces  is  calculated  during  execution  of 
the  program.  When  contact  occurs  between  an  occupant  surface  and 
a  contact  surface,  the  time  and  relative  velocity  of  impact  are 
computed  and  stored  for  printing.  The  impact  conditions  deter¬ 
mined  in  this  way  can  be  used  in  evaluation  of  injury  potential 
for  a  given  cockpit  configuration. 

4.4.2  Injury  Criteria 

The  injury  criteria  used  in  the  program  were  selected  as  the  most 
suitable  for  aircraft  crash  analysis.  SOM-hA  output  can  be  used 
to  determine  the  potential  for  injury  to  three  regions  of  the 
body;  head,  vertebral  column,  and  thorax.  Each  of  these  is  dis¬ 
cussed  in  the  remainder  of  this  chapter. 

4. 4. 2.1  Head  Injury.  An  accepted  criterion  for  head  injury  is 
the  Severity  Index  (SI)  developed  by  C.  w.  Gadd  preferences  11 
and  18),  whieh  is  calculated  for  tho  head  and  chest  according  to 

tf  „ 

SI  «  r  a  dt  (115) 

4 

where  a  “acceleration  in  G  as  a  function  of  time 
n  *  weighting  factor,  2.5  for  head  impacts 
t  =  time  in  seconds. 

Although  Gadd  used  uniaxial  acceleration  in  his  validation  of  the 
Severity  Index,  Federal  Motor  >/ehicle  Safety  Standard  208  requires 
tho  use  of  resultant  acceleration.  A  tolerable  SI  value  of  1000 
is  accepted  for  frontal  impact  of  the  head,  and  a  value  of  1500, 
for  distributed,  or  noncontact  accelerations  (reference  19). 

Also,  the  Head  Injury  Criterion  UlIC)  contained  in  Federal  Motor 
Vehicle  Safety  Standard  208  is  calculated  according  to 
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2.5 


HIC  -  max 

where  a  is  the  resultant  head  acceleration  in  G  and  and  t2  are 
any  two  points  in  tims  (sec)  during  the  crash  event,  a.  tolerable 
value  of  1000  is  accepted  for  the  HIC. 

4. 4. 2. 2  Thoracic  Injury.  Although  chest  deflection  is  commonly 
accepted  as  the  preferred  tolerance  criterion,  it  is  not  pres¬ 
ently  predicted  by  SOM -LA, 

FMVSS  208  currently  specifies  as  acceptable  an  acceleration  pulse 
whieb  "...  shall  not  exceed  60  G  except  for  intervals  whose 
cumulative  duration  is  not  more  than  3  msec,-  where  the  accelera¬ 
tion  is  the  resultant  measured  at  the  center  of  gravity.  This 
acceleration  is  printed  by  SGM-LA.  Previously,  FMVSS  208  had  ap¬ 
plied  a  Severity  Index  to  the  chest  acceleration  pulse.  This  in¬ 
dex  was  calculated  in  exactly  the  same  manner  as  the  head  Severity 
index  discussed  previously,  and  the  limit  of  1000  was  the  same  as 
that  for  the  head.  SGM-LA  also  print!*  this  ai.  tor  the  thorax, 

4. 4.2.3  Vertebral  Injury,  For  determination  of  the  potential  of 
vertebral  injury,  the  SOM-LA  two-dimensional  occupant  model  has 
beam  elements  in  the  lumbar  spine  and  neck.  Axial  forces  and 

-bending  moteents  are  output  for  these  segments,  tolerable  Corea 
levels-  for  the  lumbar  spine  hove  not  yet  been  established.  How¬ 
ever,  based  on  cadaver  tests,  references  20  and  21  have  published 
to ivrahie  bonding  moments  for  the  neck  such  as  1200  in. -lb  in 
fluxion  {forward  bending)  and  S00  iiw-Ib  ft-lb  in  extension  (rear* 
ward  -bending),  respectively. 

The  dynamic  chipoftse  index  (Dkl)  is  also  computed  by  SOM- LA  a i  a 
measure  of  the  probability  of  spinal  injury  due  to  a  vertical  ac* 
cele rat ion  parallel  to  the  spine  (reference  22). 


t  £ 
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adt 


(t2  -  t^) 


(116) 


In  this  model,  the  response  of  the  body  to  acceleration  parallel 
to  the  spine  is  modeled  by  a  single  lumped-mass,  damped-spring 
system  as  shown  in  figure  32,  or,  in  other  words,  the  total  body 
mass  that  acts  on  the  vertebral  column  to  cause  deformation  is 
represented  by  the  single  mass.  In  general,  the  motion  of  the 
system  shown  in  figure  31  obeys  the  relationship 


d2*  ^  0„  dfi  A 

+  2tU)_  TT  + 


dt 


n  dt 


n 


6  =  Z 


(117) 
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m  =  mass  (lb-sec  /in.) 
6  =  deflection  (in.) 

C  =  damping  ratio 
k  =  stiffness  (lb/in.) 

Z  =  acceleration  input 
(in. /sec2) 


*DRI 


<d_ 


n 


g 

natural  frequency  of 
the  analog  =»  \Zk/m 
(rad/sec) 


♦Dynamic  Response  Index 


g  =  386  in. /sec' 


Figure  32.  Model  used  for  prediction  of  spinal 
injury  (from  reference  22). 


The  solution,  the  deflection  6,  is  representative  of  the  deforma¬ 
tion  of  the  spine,  and  the  last  term  of  the  left-hand  side  of 
equation  (117),  divided  by  the  gravitational  acceleration,  is  the 
DRI.  The  properties  used  in  the  model  were  derived  from  tests 
involving  human  subjects  and  cadavers.  For  example,  the  spring 
stiffness  k  was  determined  from  tests  of  human  cadaver  vertebral 
segments?  damping  ratios  were  determined  from  measurements  of 
mechanical  impedance  of  human  subjects  during  vibration  and  im- 

ii 

pact.  The  acceleration  input,  Z,  is  the  component  of  seat  pan 
acceleration  parallel  to  the  vertebral  column. 
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Equation  (112)  is  solved  for  6  simultaneously  with  the  occupant 
equations  of  motion,  using  the  constant"?  x,  -  0.224  and  wn  =  52.9 
rad/sec.  The  DRI  is  then  calculated  at  each  step 

DRI  -  u)  26/g 
n  J 


(118) 


5 . 0  MODEL  VALIDATION 


Validation  has  been  based  on  data  from  several  series  of  decelera¬ 
tion  sled  tests  conducted  at  the  FAA  Civil  Aeromedical  Institute 
(CAMI).  The  response  of  the  combined  occupant  and  seat  models 
has  been  verified  by  comparison  with  data  from  tests  that  utilized 
specially  designed  and  fabricated  seats  that  had  replaceable  legs. 
Test  conditions  were  selected  to  cause  significant  plastic  deform¬ 
ation  of  the  legs.  Response  of  the  occupant  model,  particularly 
to  a  vertical  input  acceleration,  was  validated  using  data  from 
other  test  series  that  conducted  with  a  rigid  seat  and  with  a  pro¬ 
duction  energy-absorbing  helicopter  seat.  These  tests  are  de¬ 
scribee  in  further  detail  in  the  remainder  of  this  chapter. 

5.1  SIMPLE  SEAT  STRUCTURES 

Two  series  of  deceleration  sled  tests  were  performed  at  CAMI  spe¬ 
cifically  to  provide  data  for  validation  of  Program  SOM-LA.  The 
tests  utilized  an  Alderson  VIP-50  dummy  in  forward-facing  test  seats. 
The  test  program  is  described  in  reference  23,  which  includes  a 
summary  of  measured  data  from  all  the  tests. 

The  first  series  of  tests  used  a  rigid  seat  pan  and  seat  back  as¬ 
sembly,  supported  by  solid,  rectangular  cross-section  legs  and 
seat  back  hinges.  The  second  series  of  validation  tests  used  a 
similar  rigid  seat  pan  and  back,  braced  at  a  90-degree  included 
angle,  as  illustrated  in  figure  33.  The  seat  legs  were  1-in. 
diameter,  0.068-in.  wall  thickness  1010  steel  tubing,  pin  jointed 
at  the  bottom,  and  fixed  at  the  seat.  Cushions  were  1-in.  thick 
Ensolite  pads  on  the  seat  pan  and  back,  and  the  restraint  system 
consisted  of  a  conventional  nylon  lap  belt  attached  to  the  seat 
pan  with  a  double  shoulder  belt  that  was  anchored  to  the  seat 
back  and  fitted  to  the  buckle  at  the  center  of  the  lap  belt.  For 
all  tests,  the  belts  were  adjusted  to  a  snug  fit  with  all  slack 
removed.  For  each  of  these  seat  designs,  two  impact-vector 
orientations  were  used.  The  first  orientation  provided  pure 
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NOTE:  All  dimensions  in  inches. 

Figure  33.  CAMI  Series  2  test  seat. 

forward-facing  (-G  )  acceleration.  The  second  orientation  pro- 

vided  combined  longitudinal  (-G  )  and  vertical  (+G  )  acceleration 

X  z 

by  reorienting  the  seat  system  so  that  the  impact  vector  fell 
60  degrees  below  the  floor  plane  of  the  seat.  For  the  seat  de¬ 
sign  with  tubular  legs,  eight  static  tests  and  58  dynamic  tests, 
which  used  acceleration  levels  of  5.4  G  and  9.5  G,  were  conducted. 
For  the  dynamic  tests,  the  lower  acceleration  level  produced  min¬ 
imal  plastic  deformation  of  the  seat  legs  without  significant 
cross  section  change,  while  the  higher  acceleration  level  produced 
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marked  plastic  deformation  with  localized  buckling  and  cross  sec¬ 
tional  change  at  the  fixed  end.  For  the  tests  with  the  angled 
floor,  acceleration  levels  of  13.5  G  and  22  G  were  required  to 
produce  similar  results.  The  impact  velocity  for  all  of  thi^se 
tests  was  approximately  44  ft/sec.  s* 


5.1.1  Low-Deceleration  Tests 

The  first  set  of  test  conditions  considered  is  that  of  the  low- 
deceleration,  forward-facing  series.  Pre-  and  posttest  photo¬ 
graphs  from  one  of  the  10  tests  are  shown  as  figures  34  and  35. 

The  permanent  deformation  at  the  top  of  the  legs  can  be  seen  in 
the  latter  photograph.  The  trapezoidal ly  shaped  input  decelera¬ 
tion  is  shown  in  figure  36.  This  case  is  treated  in  detail  as  an 
example  in  Volume  II  -  SOM-LA  User  Manual. 

Predicted  dummy  accelerations,  restraint  system  loads,  and  forward 
displacement  of  the  seat  pan  are  compared  in  figures  37  through  41 
with  the  mean  of  data  measured  in  10  tests. 

5.1.2  Higher-Deceleration  Tests 

A  second  series  of  forward-facing  tests  used  the  sled  deceleration 
shown  in  figure  42.  As  shown  in  the  posttest  photograph  of  fig¬ 
ure  43,  a  much  greater  forward  displacement  of  the  seat  structure 
was  achieved  through  plastic  deformation  and  local  buckling  of 
the.  legs  at  their  connection  to  the  seat  pan.  Predicted  response 
is  compared  in  figures  44  through  48  with  the  mean  of  data  from 
10  tests. 

5.2  ENE RGY-ABSORBING  SEAT  TESTS 


As  part  of  a  U.S.  Army-sponsored  research  program,  a  number  of 
sled  tests  were  conducted  at  CAMI  using  a  production  crewseat  for 
the  Sikorsky  UH-60A  Black  Hawk  helicopter.  The  seat,  described 
in  detail  in  reference  24,  consists  of  a  bucket  which  supports  the 


99 


4$ 


,  g# ' 


MB  CM  0 


n4gS>U£h BM&f 

^MsSE-'i 
lEJBHiiiai 

B-lKS"1 

%& 

jHffifcfeF 

«W  S  '  f  y-'  "‘  •  "' 

UtfMs,  .:  '/-  ■ 
^  •A.v-,^. 

mm 


,  /*<«  «***  "Mil-.  , 


l_. 


figure  35.  CAM!  Series  2  „ 

iorwerd-facln' 


101 


Forward  deflection,  in 


Time,  sec 


Figure  41.  CAMI  validation  Series  2,  low-deceleration,  forward 
facing  tests,  seat  forward  displacement. 


TiB«  -  »»ec 


Figure  42.  Sled  deceleration,  CAMI  Series  2,  higher- 
deceleration,  forward-facing  tests. 
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(a)  Lap  belt 


(b)  Shoulder  belts 

Figure  47.  CAMI  validation  Series  2,  higher-deceleration,  forward¬ 
facing  tests,  restraint  system  loads. 
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Figure  48.  CAMI  validation  Series  2,  higher-deceleration,  forward¬ 
facing  tests,  seat  forward  displacement. 

occupant  and  to  which  the  restraint  system  is  attached  and  a  frame 
mounted  on  the  aircraft  floor.  Principal  functional  members  of 
the  frame  are  two  vertical  (or  nearly  vertical)  guide  tubes  along 
which  the  bucket  can  move,  controlled  by  one  or  more  energy¬ 
absorbing  devices.  Vertical  inertial  crash  loads  force  the  seat 
bucket  down  the  guide  tubes  against  the  resistance  of  the  energy 
absorbers,  producing  an  energy-absorbing  stroke  in  that  direc¬ 
tion.  For  most  efficient  use  of  the  stroke  distance  available 
between  the  bucket  and  the  floor,  energy  absorbers  are  designed 
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to  stroke  at  constant  load,  that  load  being  determined  by 
design  criteria  based  on  human  tolerance  to  +G^  acceleration. 

In  order  to  simulate  a  nearly  vertical  impact  on  the  horizontal 
sled,  the  seat  z-axis  was  pitched  forward  17  degrees  from  the 
plane  of  the  sled,  as  shown  in  figure  49.  A  pretest  view  of  the 
dummy  positioned  in  the  seat  is  presented  in  figure  50,  where  the 
fixture  simulating  the  aircraft  floor  can  be  seen.  (The  Black 
Hawk  helicopter  has  a  well  beneath  each  crewseat  to  permit  addi¬ 
tional  stroke,  so  that  a  minimum  of  12  in.  and  a  maximum  of  17  in. 
can  be  attained,  depending  on  the  seat's  vertical  adjustment  posi¬ 
tion.  )  The  dummy  was  flexed  forward  at  the  waist  for  positioning 
in  the  seat,  and  its  lower  torso  was  pushed  firmly  against  the 
back  cushion.  The  five-point  restraint  system  was  installed,  and 
the  lap  belt  was  tightened.  With  the  inertia  reel  unlocked,  the 
shoulder  harness  was  tightened  using  its  adjusters,  and  a  preload 
was  applied  to  the  reel.  The  dummy's  feet  were  taped  to  the  foot¬ 
rest,  which  was  supported  on  a  six-axis  load  cell.  For  the  test 
used  in  SOM-LA  validation,  the  impact  velocity  was  43.5  ft/sec, 
and  the  sled  deceleration,  with  a  peak  of  41.5  G,  was  as  shown  in 
figure  51. 

The  pelvic  structure  of  the  50th-percentile  Part  572  dummy  used 
in  the  test  had  been  modified  to  include  a  six-axis  load  cell  at 
the  base  of  the  lumbar  spine  (reference  25).  Forces  were  thus 
measured  in  the  dummy  spine  as  well  as  in  the  footrest.  Acceler¬ 
ations  were  measured  in  the  dummy  pelvis,  chest,  and  head,  and 
on  the  seat.  Seat  stroke  was  also  measured  using  a  displacement 
transducer. 

The  two-degree-of-freedom  seat  model  was  used  in  simulating  this 
test.  The  energy  absorber  force-deflection  characteristics  were 
based  on  data  from  static  and  dynamic  tests  of  these  components. 
The  seat's  rotational  stiffness  was  based  on  a  static  test  of  the 
system  (reference  24). 
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Figure  49.  Test  configuration  for  CAMI  tests  with 
energy-absorbing  helicopter  seat. 

Comparisons  of  predicted  accelerations,  seat  stroke,  and  forces 
with  measured  data  are  presented  in  figures  52  through  61. 

5.3  RIGID  SEAT  TESTS 


As  part  of  the  same  U.S.  Army-sponsored  program  mentioned  in  sec¬ 
tion  5.2,  a  number  of  tests  were  conducted  with  various  dummies  in 
a  rigid  seat  whose  seat  pan  and  back  formed  a  right  angle  with  re¬ 
spect  to  each  other,  as  shown  in  figure  62  (where  a  95th-percentile 
dummy  is  shown).  No  cushions  were  used,  and  the  plywood  seat  pan 
was  supported  by  a  six-axis  load  cell.  The  four-point  restraint 
system  used  automotive -type  nylon  webbing,  the  vertical  (z)  axis 
of  the  seat  was  aligned  with  the  velocity  vector,  in  other  words, 
horizontal  for  the  sled  impact.  For  the  test  with  the  modified 
Part  572  dummy,  the  sled  deceleration  is  shown  in  figure  63. 


0  0.05  O.in  0.15  0.20  0.25 

Time,  sec 

Figure  51.  Sled  deceleration  for  energy-absorbing  seat  test. 

Accelerations  were  measured  in  the  dummy,  and  forces  were  measured 
in  the  lumbar  spine,  footrest,  and  seat  pan. 

Simulation  results  and  test  data  are  compared  in  figures  64  through 
70. 

5.4  GENERAL  AVIATION  SEAT  TESTS 

To  complete  the  SOM-LA  validation,  CAMI  conducted  several  teats 
with  production  general  aviation  seats,  which  were  purchased  from 
dealers  of  used  aircraft  hardware.  Twelve  seats  of  two  different 
designs  were  used  in  a  total  of  four  static  and  thirteen  dynamic 
tests.  The  tests  which  will  be  described  here  utilized  a  seat 
with  a  frame  of  rectangular  aluminum  tubing.  This  seat  produced 
more  interesting  results  for  validation  of  the  seat  model  than  the 
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Figure  52.  Energy-absorbing  helicopter  seat  test, 
seat  z-acceleration. 
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Figure  54.  Energy-absorbing  helicopter  seat  test 
dummy  chest  z-acceleration. 


CHEST  Z-fiXIS  RCCELERfiTION  C Q) 


100 


-  Simulation 

-  Test  data 


Figure  54.  Energy-absotbing  helicopter  seat  test 
dummy  chest  z-acceleration. 
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Figure  SS.  Energy-absorbing  helicopter  seat  test 
duraoy  head  z-acceleration. 


Figure  56.  Energy-absorbing  helicopter  seat  test, 
seat  vertical  displacement. 
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Figure  58.  Energy-absorbing  helicopter  seat  test, 
tootrest  x-force. 


FOOTREST  Z-RXIS  FORCE  (LB) 


Figure  59.  Energy-absorbing  helicopter  seat  test 
footrest  z-force. 
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V” 


TIME  (SEC) 

Figure  60.  Energy-absorbing  helicopter  seat  test, 
lumbar  spine  axial  load. 


LUMBRR  SPINE  MOMENT  (IN-LB) 


Figure  61.  Energy-absorbing  helicopter  seat  test, 
lumbar  spine  moment. 


127 


82  09003  61 


Figure  62.  CAMI  rigid  seat  (with  VIP-95  dummy  in  place). 
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Figure  64. 


Rigid  seat  test,  dummy  pelvis  z-acceleration 
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Figure  66.  Rigid  seat  test,  dummy  head  z-acceleration. 
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82  09003  68 


CUSHION  FORCE  (LB) 


Figure  67.  Rigid  seat  test,  seat  pan  z-force. 
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Figure  69.  Rigid  seat  test,  lumbar  spine  axial  force. 
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second  seat  because  of  bending  in  the  structure.  The  seat  and 
50th-peccentile  (Part  572)  dummy  prior  to  the  initial  dynamic  test 
ore  shown  in  figure  71. 

In  order  to  achieve  sufficiently  high  decelerative  loading  to  cause 
deformation  of  the  seat  frame,  it  was  necessary  to  modify  the  floor 
(sled)  attachments.  The  modified  right-front  leg  of  the  seat  is 
shown  in  figure  72,  where  it  can  be  seen  that  the  track  roller 
has  been  removed  and  the  track  roller  housing  has  been  pinned  to 
an  inverted  steel  T-section.  The  web  of  the  T-section  has  been 
notched  to  permit  rotation  of  the  seat  leg.  The  two  rivets  that 
secure  the  roller  housing  to  the  lower  end  of  the  leg  have  been 
replaced  by  0.25-in.  steel  bolts.  All  four  of  the  seat  legs  were 
modified  in  this  manner. 

Each  leg  was  bolted  to  a  six-axis  load  cell.  A  general  aviation- 
type  three-point  restraint  system  was  used,  with  the  shoulder  belt 
passing  over  the  dummy's  left  shoulder  and  anchored  above  and  to 
the  left  of  the  seat,  in  a  position  typical  of  an  aircraft  installa¬ 
tion.  Forces  were  measured  in  both  sides  of  the  lap  belt  and  in 
the  shoulder  belt,  between  the  dummy  and  the  anchor  point.  As 
shown  in  figure  71,  a  displacement  transducer  was  installed  to 
measure  forward  motion  of  the  seat  frame.  For  some  tests  a  foot¬ 
rest  was  included,  mounted  at  an  angle  of  45  degrees  relative  to 
the  horizontal  and  bolted  to  a  load  cell.  In  tests  that  used  a 
flat  floor,  the  foot  forces  were  not  measured,  as  the  feet  leave 
the  floor  a  very  short  time  after  impact.  Accelerations  were 
measured  in  the  pelvis,  chest,  and  head  of  the  dummy.  The  dummy 
had  been  modified  to  include  a  load  cell  at  the  base  of  the  lumbar 
spine,  as  in  the  tests  described  in  sections  5.2  and  5.3. 

The  initial  finite  element  model  of  the  seat  structure  illustrated 
in  figure  73  used  28  nodes,  36  beam  elements,  and  2  triangular  plate 
elements.  The  approximation  to  the  cross  section  of  the  6061-T6 
aluminum  beam  elements  is  shown  in  figure  74. 


Figure  71.  General  aviation  seat  and  Part  572  dummy 
prior  to  CAMI  Test  A81-108. 


The  tests  <jiven  greatest  attention  in  the  simulations  because  of 
bending  in  the  seat  structure  were  those  conducted  in  a  forward¬ 
facing  configuration  with  a  12-G  sled  deceleration  (figure  75). 

Two  tests,  A81-110  and  A81-111,  were  conducted  under  similar  con¬ 
ditions  except  that  the  former  included  a  footrest,  as  shown  in 
figure  76,  and  the  latter  did  not  (figure  77).  Results  of  both 
tests  were  similar  except  that  without  the  footrest  dummy 
x-accelerations  reached  somewhat  higher  magnitudes.  In  test 
A81-111,  without  the  footrest,  the  shoulder  belt  load  reached  a 
higher  value,  and  the  dummy,  rebounding  rearward,  caused  failure 
of  the  seat  back.  Because  of  the  latching  mechanism  on  the  fold¬ 
ing  seat  back,  the  actual  seat  was  less  able  to  resist  rearward 
loading*  than  the  finite  element  model,  for  which  the  frame  was 
assumed  to  be  continuous  from  the  seat  pan  to  the  seat  back.  Ex¬ 
amining  the  model  in  figure  73  and  using  the  right-hand  side  of  the 
seat  as  an  example,  the  latch  is  represented  by  beam  elements  con¬ 
necting  nodes  10  to  11  and  II  to  12.  The  cross  section  of  these 

elements  is  that  of  the  tubing  shown  in  figure  74,  rather  than  the 
thinner  material  actually  used  in  the  latch.  This  structure  is 
further  supported  by  a  triangular  plate  connecting  nodes  8,  10, 
and  12,  resulting  in  a  structure  that  is  much  stronger  than  the 
actual  hardware.  In  the  SOM-LA  simulations,  the  impact  of  the 
rebounding  dummy  into  the  relatively  rigid  seat  back  produces  a 
spike  in  the  segment  accelerations,  particularly  the  chest,  as 
shown  in  figures  78  and  79.  The  effect  of  this  impact  with  the 
seat  back  on  occupant  motion  can  be  seen  in  the  last  two  plots  of 
figure  80. 

The  revised  finite  element  model  is  shown  in  figure  81.  The  short 

beam  elements  that  connect  nodes  8  to  10  and  22  to  24  contain  tor¬ 

sional  releases  to  allow  pivoting  of  the  seat  back  about  the  lateral 
horizontal  line  from  node  8  to  node  22.  Crushing  or  buckling  of 

*NAS  809,  paragraph  4. 1.2.1,  requires  the  capability  of  the  seat  to 
resist  a  300- lb  load  applied  aftward  8  in.  above  the  intersection 
of  the  seat  back  and  bottom. 
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Sled  deceleration  (G) 


Figure  75.  Sled  deceleration,  CAMI  general 
aviation  seat  Test  A81-110. 
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CHEST  X-RXIS  RCCELERRTION  (G) 


Figure  78.  General  aviation  seat  Test  A81-110, 
dummy  chest  x~accele ration  (initial 
finite  element  model). 
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Simulation 


0.0$  0.10  0.15  0.20  0.2$ 


TIME  (SEC) 


Figure  79.  General  aviation  seat  Test  A81-111, 
dummy  chest  ^-acceleration  (initial 
finite  element  model)  > 


PROGRAM  SOM-LA  OCCUPANT  MOOEL 
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TIME-  .1600  SEC. 


PROGRAM  SOM-LA  OCCUPANT  MOOEL 
NONAOJUSTA6LE  PILOT  SEAT  12-0  FORWARD  TEST  A81-110 
TIME  -  .2000  SEC. 


PROGRAM  SOM- LA  OCCUPANT  MODEL 
KONAOJUSTflHLE  PILOT  SEAT  12-0  FORWARD  TEST  A61-UQ 
TIME  -  .2400  SEE. 


PSORA  AM  SOM-LA  OCCUPANT  MCCCL 
NQKADJUSTABLC  PILOT  SEAT  13-0  rORWARO  TEST  ADl-lIO 
TIME  -  .2800  SEC. 


■n- 


M\ 


Figure  80  (contd).  General  aviation  seat  Test  A81-110, 

predicted  occupant  position  {initial 
finite  element  model). 
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PROGRAM  SOM-Lfl  SEAT  STRUCTURE  MODEL 
NONADJUSTRBLE  PILOT  SEAT  12-8  FORWARD  TEST  A81-110 
PLOT  NO.  1,  TIME  -  0.0000  SEC. 


Figure  81.  Revised  finite  element  model  of 
general  aviation  seat. 


the  seat  back  latch  mechanisms  is  modeled  by  nonlinear  spring  ele¬ 
ments  that  connect  nodes  11  to  12  and  25  to  26.  A  complete  list¬ 
ing  of  the  input  data  is  presented  in  figure  82. 

The  revisions  to  the  seat  model  reduced  the  severity  of  the  occupant' 
impact  with  the  seat  back  by  permitting  deformation  of  the  latches 
(spring  elements).  However,  at  0.237  sec,  after  nodes  9  and  23 
at  the  top  of  the  seat  back  had  been  pushed  more  than  5.4  in.  rear¬ 
ward,  the  simulation  was  halted  due  to  ill  conditioning  in  the 
seat  stiffness  matrix.  Occupant  position  just  prior  to  this  time 
can  be  seen  in  figure  83,  and  the  seat,  in  figure  84. 

The  remaining  SOM-LA  predictions  for  test  A81-110  (45-degree  foot¬ 
rest)  are  compared  with  measured  data  in  figures  85  through  100. 
Footrest  forces  are  included  (figures  97  and  98),  although  the  X- 
and  Z-  directions  for  simulation  and  test  do  not  coincide.  SOM-LA 
assumes  the  X-direction  in  the  plane  of  the  floor  whereas  for  these 
tests  the  load  cell  was  rotated  45  degrees  (figure  71).  A  transfor¬ 
mation  between  the  coordinate  systems  would  thus  be  required  for 
comparison  of  forces.  Lumbar  axial  force  and  y-moment  measured  in 
the  test  are  included  (figure  99  and  100).  However,  because  of  the 
nonsymmetric  three-point  restraint  system,  the  three-dimensional 
occupant  model  was  used,  and  this  model  does  not  have  the  capability 
for  spinal  load  prediction. 

5.5  OISCUSS ION 

The  simulation  results  presented  in  this  chapter  were  obtained 
using  identical  modeling  parameters  in  all  cases.  For  any  single 
case,  performance  could  be  improved  by  adjustment  of  certain  in¬ 
put  parameters,  but  the  standard  values  as  presented  in  this  re¬ 
port  (including  Volume  II)  are  based  on  compromises,  c.g.,  between 
vertical  and  longitudinal  impact  configurations. 

The  most  likely  parameter  for  optimizing  seat  response  in  a  parti¬ 
cular  case  is  the  buckling  coefficient,  described  in  section  3.6. 
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Figure  82.  Input  data  listing  for  simulation 
of  Test  A81-110. 
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Figure  82  (contdl.  Input  data  listing  for  simulation 

of  Test  A81-U0. 
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PROGRAM  SOM-Lfl  OCCUPANT  MODEL 
NONADJUSTRBLE  PILOT  SEAT  12-8  FORWARD  TEST  A81-U0 
TIME  -  .2350  SEC. 
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Figure  83.  General  aviation  seat  Test  A81-110, 
predicted  occupant  position  at 
t  -  0.235  sec  (revised  finite 
element  model). 
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PROQRRM  SOM-Lfi  SERT  STRUCTURE  MODEL 
NQNRDJUSTRBLE  PILOT  SERT  12-0  FORWRRD  TEST  R81-110 
PLOT  NO.  8,  TIME  -  .2350  SEC. 
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Figure  84.  General  aviation  seat  Test  A81-110,  predicted 
seat  position  at  t  =  0.235  sec  (revised  finite 
element  model). 


Figaro  85.  General  aviation  seat  Test  A31-110, 
pelvis  x-acceleration. 
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Figure  86.  General  aviation  seat  Test  A31-110, 
pelvis  z-acceleration. 
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Figure  87.  General  aviation  seat  Test  A81-110 
pelvis  resultant  acceleration. 


CHEST  X-RXIS  ACCELERATION  (G) 


Figure  88.  General  aviation  seat  Test  A81-110 
chest  x-acceleration. 
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Figure  89.  General  aviation  seat  Test  A81-110 
chest  s-acceleration. 
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Figure  91.  General  aviation  seat  Test  A81-110, 
head  x-acceleration. 
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Figure  92.  General  aviat 
head  z-accele 
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Figure  93.  General  aviation  seat  Test  A81-110, 
head  resultant  acceleration. 
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Figure  95.  General 
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Figure  99.  General  aviation  seat  Test  A81-110, 
lumbar  axial  force. 
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Figure  100.  General  aviation  seat  Test  A81-110 
lumbar  y-moment. 


The  value  of  0.5,  which  has  a  rational  basis,  produced  seat  dis¬ 
placement  in  the  low-deceleration  case  which  was  within  the  en¬ 
velope  of  the  test  data  (figure  41).  For  the  higher-deceleration 
case,  the  predicted  seat  displacement  was  low  (figure  48),  but 
could  be  increased  by  increasing  the  coefficient.  However,  another 
reason  for  the  reduced  seat  displacement  could  be  the  model's  ig¬ 
noring  seat  inertia.  In  this  case  the  seat  weight  was  more  than 
50  lb,  so  that  its  contribution  to  the  displacement  measured  in 
the  test  may  have  been  significant.  However,  general  aviation 
seats  are  expected  to  weigh  considerably  less,  so  that  this  dis¬ 
crepancy  should  be  reduced  in  simulation  of  actual  seat  systems. 

Occupant  parameters  which  might  be  adjusted  to  improve  simulation 
of  one  particular  type  of  test  configuration,  such  as  longitudinal 
impacts,  are  the  chest  and  abdominal  compliance  parameters,  which 
are  used  to  effectively  soften  the  input  belt  characteristics. 

Also,  if  head  response  is  of  particular  interest,  the  neck  damp¬ 
ing  coefficients  can  be  adjusted  to  achieve  optimum  response. 
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6.0  CONCLUSIONS 


The  Seat/Occupant  Model  -  Light  Aircraft  (Program  SOM-LA)  provides 
a  method  for  analyzing  the  performance  of  a  seat  and  occupant  in 
an  aircraft  crash  environment.  For  preliminary  analyses  of  cabin 
configuration  or  restraint  system  design  where  details  of  the  seat 
structure  are  not  known,  a  rigid  seat  can  be  used.  Otherwise, 
finite  element  modeling  capability  is  available  for  the  seat 
structure. 

Modeling  parameters  have  been  based  on  a  compromise  among  various 
impact  configurations  and  seat  types.  If  a  user  were  to  desire 
more  rigorous  modeling  of  a  particular  impact  configuration,  cer¬ 
tain  parameters  could  be  adjusted  to  optimize  simulation. 
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APPENDIX  A 

OCCUPANT  SEGMENT  POSITION 
THREE-DIMENSIONAL  MODEL 


Referring  to  figure  A— 1,  and  using  p^  =  —  p^  the  absolute 

position  of  the  mass  center  of  each  body  segment  is  given  below. 
The  elements  of  the  transformation  matrices  [Tn]  are  functions 
of  the  generalized  coordinates,  as  given  by  equations  (6),  (9), 
and  ( 10) . 

Segment  1: 

(Xr  Z^),  the  coordinates  of  the  reference  point  on  the  body 
are  the  generalized  coordinates  (q^  q2,  q3). 


Segment  2: 


Segment  5: 


Segment  6: 


Segment  7 

1 


Segment 


t 

appendix  b 

OCCUPANT  SEGMENT  POSITION; 
TWO-DIMENSIONAL  MODEL 


Referring  to  figure  6,  the  absolute  position  of  the  mass  center 
of  each  body  segment  is  given  below.  Lg  and  LR  are  the  lateral 
distances  from  the  mid-saggital  plane  to  the  shoulder  and  hip 
joints,  respectively. 


Segment  1: 

(X  Y  ,  Zt) ,  the  coordinates  of  the  reference  point  on  the  body. 
X1  and  are  the  generalized  coordinates  q x  and  q2. 


Segment  2: 


S  (cos  0.  -  cos  02) 
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a2  3  +  el  sin  °1 
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Segment  3; 
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Segment  4: 


X4  =  Xi  -  el  cos  01 
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Segment  5: 


X5  =  Xt  -  eL  cos  Ql 
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Segment  6: 
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